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16. Abstract 

Fixed- time-of -arrival (FTA) guidance and navigation is investigated as a possible technique 
capable of operation within much more stringent en route separation standards and offering signifi 
cant advantages in safety, higher traffic densities, and improved scheduling reliability, both en 
route and in the terminal areas. 

This study investigates the application of FTA guidance previously used in spacecraft 
guidance. These FTA guidance techniques have been modified and are employed to compute the 
velocity corrections necessary to return an aircraft to a specified great-circle reference 
path in order to exercise en route time and position control throughout the entire flight. 

The necessary position and velocity estimates to accomplish this task are provided by Kalman 
filtering of data from Loran-C, VORTAC/TACAN, Doppler radar, radio or barometric altitude, 
and altitude rate. The guidance and navigation system was evaluated usi^ig a digital simula- 
tion of the cruise phase of supersonic and subsonic flights between San Francisco and New 
York City, and between New York City and London. 

It was found that the £ixed-time-of- arrival guidance system is capable of easily meeting 
the present lateral separation standards and of maintaining en route timing errors to less 
than 30 sec over intercontinental routes. 
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SYMBOLS 


d slant range from transmitter to receiver 

H Jacobian matrix relating data-set quantities and vehicle state 

(2x6 matrix) 

I identity matrix with appropriate dimensions 
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~ ~T 
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E(xx^) 

AP correction to the P matrix following a course correction maneuver 

As commanded change in ground speed 

S vector from center of earth to a station on surface 

SNR signal-to-noise-ratio 

t general time argument, hr 

time of the kth data-set sample 

t7 time of the kth data-set sample but further updating required 

K 

St deviation of the arrival time from the reference arrival time 

v(t) velocity deviation vector from its reference value 
V total velocity vector, with respect to earth-fixed coordinate frame 
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EN ROUTE POSITION AND TIME CONTROL OF AIRCRAFT USING 
KALMAN FILTERING OF RADIO AID DATA 
Leonard A. McGee and Jay V. Christensen 
Ames Research Center 


SUMMARY 


Fixed-time-o£-arrival CPTA) guidance and navigation is investigated as a 
possible technique capable of operation within much more stringent en route 
separation standards and offering significant advantages in safety, higher 
traffic densities, and improved scheduling reliability, both en route and in 
the terminal areas. 

Thus study investigates the application of FTA guidance previously used 
in spacecraft guidance. These FTA guidance techniques have been modified and 
are employed to compute the velocity corrections necessary to return an air- 
craft to a specified great-circle reference path in order to exercise en route 
time and position control throughout the entire flight. The necessary posi- 
tion and velocity estimates to accomplish this task are provided by Kalman 
filtering of data from Loran-C, VORTAC/TACAN , Doppler radar, radio or baro- 
metric altitude, and altitude rate. The guidance and navigation system was 
evaluated using a digital simulation of the cruise phase of supersonic and 
subsonic flights between San Francisco and New York City, and between New York 
City and London. 

Performance is evaluated with emphasis on the ability of the system to 
maintain position and time control along the great-circle reference paths. 

It was found that the fixed-time-of-arrival guidance system is capable of 
easily meeting the present lateral separation standards and of maintaining 
en route timing errors to less than 30 sec over intercontinental routes. 


INTRODUCTION 


As air traffic density steadily increases, it is important to improve 
accuracy in the estimation and control of aircraft position and velocity so 
that separation of aircraft may be reduced without increasing the risk of col 
lision. In particular, if a reliable and accurate en route timing capability 
can be attained, then significant advantages in safety, higher traffic densi- 
ties, air traffic control procedures, and economies associated with aircraft 
scheduling will result. 

It is felt that on-board digital computer technology has improved to the 
point where advanced computer-oriented integrated systems techniques can be 



considered in the improvement of aircraft navigation and guidance. The 
larger, more expensive jet aircraft, such as the jumbo jets and supersonic 
transports may well find the application of a computer-oriented system to be 
economically justifiable. Such a system would hold potential for improvement 
of the present method of arrival time control. 

This study investigates the performance of a computer-oriented system 
which utilizes Kalman filtering of data from VORTAC/TACAN , Loran-C, a radio 
altimeter, a barometric altimeter, and Doppler radar for navigation. The 
guidance technique of this investigation is based on the well-known fixed-time- 
of-arrival guidance previously used successfully for midcourse guidance of 
spacecraft (ref. 1) . This guidance technique is modified appropriately for 
the aircraft cruise phase guidance problem and is used by the navigation and 
guidance system throughout the flight to compute the necessary changes in 
direction and speed that will return the aircraft to a reference great-circle 
route and maintain the aircraft on-time throughout the entire flight. The 
results obtained are from a simulation on a digital computer of the cruise 
phase of flight for supersonic and subsonic flights between San Francisco and 
New York City, and between New York City and London. The necessary position 
and velocity estimates are provided by Kalman filtering of data from Loran-C, 
VORTAC/TACAN, Doppler radar, a radio altimeter, or a barometric altimeter, 
trror models for these subsystems are derived and discussed. The radio aids 
were assumed to be of the latest types with ground equipment locations cur- 
rently operational along the selected routes. The following factors were 
included in the study: (a) attitude measurement uncertainty; (b) aircraft 

wander between data samples due to external disturbances; and (c) number and 
timing of course corrections. 


Performance is evaluated with emphasis on the ability of the system to 
meet an established time of arrival at a fixed position. Performance is also 
evaluated on the basis of: (a) error in the estimation of position and veloc- 

ity; (b) course correction requirements in terms of heading and velocity mag- 
nitude changes that are required; and (c) the ability of the system to control 
the aircraft position and velocity with respect to a selected great-circle 
reference route. 


SYSTEM DESCRIPTION 


General Description of the On-Board System 

The on-board navigation and guidance system that was investigated is 
illustrated in figure 1. Simulated data from the sensors included random 
errors generated from sensor error models. 

The functions carried out by the estimation and guidance system computer 
are receiver channel selection and sensor sampling, data filtering and esti- 
mation of the aircraft state, and course correction commands for fixed- time- 
o -arrival guidance. The on-board sensors sampled successively by the com- 
puter for data filtering are the Doppler radar system, which measures 
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Attitude and occcierqlion 



Figure 1.— The on-board navigation and guidance system. 

groundspeed and drift angle, and the radio altimeter or the air-data system, 
both of which measure altitude and altitude rate. The computer also samples 
other sensors used in data filtering computations which derive their basic 
information from ground station radio transmissions. The Loran-C receiver is 
used to obtain signals from which latitude and longitude are computed by a 
coordinate converter. The VORTAC/TACAN receivers are used to determine the 
range and bearing of the aircraft relative to the selected ground station. 
Other on-board sensors which are not used directly by the computer in filter 
computations but are used by the other on-board systems are the vertical and 
heading references, which measure aircraft pitch, roll, and heading, and a 
body-mounted accelerometer, which measures longitudinal acceleration. 

Also on-board is an automatic control system which is in operation con- 
tinuously throughout all flights. Its functions are to: (1) change the air- 

craft heading in accordance with computer commands necessary to fly a great- 
circle course; (2) maintain heading within specified limits; (3) maintain 
altitude within specified limits; and (4) respond to course correction com- 
mands from the on-board computer. This system cannot fly the aircraft along 
a preprogrammed reference path without error because the on-board sensors on 
which it must rely for information are imperfect measuring devices. As a 
result of these imperfections and external disturbances, the aircraft will 
drift from the reference course to such an extent that some corrective action 
is necessary for both en route safety and adherence to schedule. This cor- 
rection is provided by an on-board navigation system, which determines the 
aircraft position and velocity, and the on-board guidance system, which com- 
putes the necessary corrections to the groundspeed, altitude, and heading 
which will return the aircraft to the reference path. Indeed, a sufficient 
number of course corrections are desirable to keep the aircraft within a 
specified distance of the reference path, to keep it on-time, and to maintain 
safe separation from other aircraft. 
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Flight Environment 


Gveat-Q-lrole reference routes— The great-circle reference routes between 
San Francisco (SF) and New York City (NYC), and between New York City and 
London (L), are illustrated in figure 2. This figure also shows the existing 
Atlantic area Loran-C stations that were used in this study. 

The principal difference between these two flight paths is the type of 
position data that is available. For the first 75 percent of the land route 
from San Francisco to New York City, only VORTAC/TACAN is available. During 
the remaining 25 percent, both VORTAC/TACAN and Loran-C are available. Posi- 
tion information in the vertical direction is obtained from an on-board baro- 
metric altimeter which is part of an air-data system and is used in prefer- 
ence to a radio altimeter, which measures the altitude above the local terrain 
instead of above the mean sea level. Over the water route from New York City 
to London, longitude and latitude information is obtained from Loran-C exclu- 
sively, even when near New York City, because it provides more accurate infor- 
mation than VORTAC/TACAN. Altitude information is obtained from a radio 
altimeter since there are no terrain effects to consider. 

Cru/vse oond/vtioTLs The cruise phases of flight for two different air- 
craft are considered. Each aircraft attempts to follow a reference great- 
circle path from origin to destination, with a fixed-time-of-arrival at the 
destination. The entire reference cruise flight from origin to destination is 
at the specified cruise altitude and groundspeed for each aircraft. The two 
reference cruise conditions considered will be referred to hereafter as 
supersonic" and "subsonic." The supersonic cruise conditions represent an 
aircraft flying at a reference groundspeed of 3,000 km/hr (1,620 knots) at 
21.336 km (70,000 ft) altitude. For the subsonic cruise, the reference ground- 
speed is 1,000 km/hr (540 knots) and the altitude is 10.668 km (35,000 ft). 


The initial and final position coordinates of the great-circle reference 
routes are summarized below: 


Latitude 

Longitude 


San Francisco 

37.634“ north 
122.359“ west 


New York 

40.634“ north 
73.783“ west 


London 

52.000“ north 


1 


.000“ west 


External disturbances— The forecast wind over the originating airport was 
assumed arbitrarily to be 100 km/hr from 210“ relative to true north. The 
random component was realized from a normal distribution with zero mean and 

along-track and cross -track directions 
ot 25 km/hr. The initial heading of the aircraft was altered to "crab" the 
aircraft into the average wind in order to aline the ground velocity vector 
along the desired track; the estimated groundspeed was adjusted accordingly. 

In this study the aircraft was considered to have sufficient reserve power so 
that there would be no difficulty in adjusting the airspeed to obtain the 
desired groundspeed. During the remainder of the flight, the effects of 
random winds were brought into the problem by random time-correlated varia- 
tions in the aircraft ground track, heading and velocity. 





Changes in the steady-state wind components were found by averaging the 
Doppler drift angle and integrating the acceleration from the body-mounted 
accelerometer along the aircraft longitudinal axis. The Doppler drift angle 
measures the angle between the ground track velocity vector and the aircraft 
longitudinal axis in the horizontal plane. Averaging this angle over several 
seconds gives an estimate of the change in the cross-track steady-state wind; 
this value was used by the computer to compensate for the change in aircraft 
heading. Similarly, integrating the output of the accelerometer (compensated 
for deviations from the level state) gives an estimate in the along-track 
change in the steady-state wind. 


Desired time en route— The main objective of this study was to investi- 
gate fixed-time-of-arrival guidance and navigation. The objective of the 
guidance system was to fly a great-circle reference course with the en route 
times given in the following table. 


Route 

Flight 

time (min) 

Subsonic 

Supersonic 

S.F. to N.Y.C. 
N.Y.C. to London 

250.079 

330.066 

83.360 

110.022 


For simplicity, a fixed sample rate was used, with adjustments to the 
rate when useful. 


Dai^ samgl'vng- For the supersonic case, sensors are selected and sampled 
every minute except for one special case, in which the sample rate is 
increased to 2/min just prior to course corrections. For the subsonic case 
sensors are sampled every 3 rain. The sensors are sampled in a sequential pat- 
tern. A typical pattern is VORTAC/TACAN or Loran-C first; Doppler radar 
second; and an altimeter, third. 


The following logic was used for accepting VORTAC/TACAN and Loran-C data- 
set s^ples: slant ranges to all the stations are computed based on the air- 

^ position and the coordinates of the stations. For VORTAC/ 

TACAN, the station with the minimum slant range is used if the slant range is 
less than 275 km (148 nmi) . If the slant range is greater than 275 km, the 
s^ple is skipped. For Loran-C, the signals from three stations are required. 

having the smallest slant range to its farthest station is selected 

accepted if: (1) the slant range is less than 

jreatS line-of-position (lOP) crossing angle is 

greater than 25 . If the slant range is greater than 2800 km, then no other 
station triad can meet the acceptance requirements and the sample is skipped. 

^ angle is less than 25®, the acceptance criteria are 

repea e or another station triad with the next larger slant range to its 

process is repeated until a suitable station triad is 
or e sample is skipped because the slant range requirement cannot be 

iner • 
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Coui'se oovpsot'ions — In selecting a course correction strategy, it 
required that the strategy demonstrate effectiveness both in reducing t e at 
eral position dispersion during flight (air corridor) and in maintaiwng t e 
aircraft on-time with a reasonable amount of aircraft maneuvering. The strat- 
egy should also function satisfactorily with thresholds which must be excee e 
by the commanded changes before execution would be allowed. These thresholds 
were imposed to reduce the probability of making course corrections whic 
would cause greater course error than already present due to errors in execut- 
ing the course correction. For the selected strategy each course correction 
consists of two maneuvers as illustrated in figure 3. The first maneuver at 



Figure 3.— Typical two-maneuver course correction. 

t\ is designed to cause the aircraft to intercept the reference flight path at 
point I after a fixed time At. The value of At is taken as 5 and 15 min, 
respectively, for the supersonic and subsonic cases. Because of external dis- 
turbances during At and errors in making the correction at t\, the aircraft 
will not exactly reach the intercept point I at t 2 - This situation is to be 
expected and, if the state estimate at t£ is reasonably accurate, causes no 
difficulty. A second correction maneuver is made at tz to intercept the des- 
tination at the desired time, tj. Again, because of control errors, naviga- 
tion errors, and course correction errors, the actual path will not exactly 
follow the predicted path as is shown in figure 3. Even though these cor- 
rection maneuvers are executed with some error, it will be shown later that 
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this is an effective course correction strategy capable of meeting the 
requirements specified. 

Course corrections are implemented as a change in aircraft speed and 
heading. A change in altitude rate is also computed but the changes were 
always less than the threshold value and therefore not used. 

To gain insight into the fundamental properties of the two-maneuver 
course correction strategy, a simplified example involving only speed and 
heading change will be examined. Consider sketch (a), which depicts the first 
maneuver of the two-maneuver strategy. 



where 

error 
track error 
reference groundspeed 
estimated groundspeed 
estimated position 

Af required heading change from V„: AH'] + A'l'o 

As required change in groundspeed 

At fixed time to intercept point 

From the geometry of sketch (a) : 

(1 7^1 + As)At = [(|7^|At + AZ^y)2 + 

= I 7^ I At [(1 + AX^y/|7^|At)2 + (A;f^y/[7^|At)2]^/2 


CD 
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from which 


tan Aipi = Aij/i 



As = ([7^1 - |7^1) + AX^y/A* 


C4) 

(5) 


At = At2 + 


( 6 ) 


From equations (S) and (6), it is seen that the speed and heading changes 
each consist o£ two terms; the first term corrects for the current error in 
velocity or heading ~ I^S”! A 1 P 2 ) second term corrects for cur- 

rent position errors. Considering position errors to a first order approxi- 
mation, the speed change. As, corrects the along-track (or en route timing) 
error and the heading change corrects the cross-track error. These are impor- 
tant relationships which will aid in interpreting the course correction data 
to be presented in a following section. 

The times at which the course corrections are scheduled for the various 
cases are given in table 1. A fixed course correction sequence is used for 
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each case to aid in comparing simulation results. Three course corrections 
are used on all the supersonic flights. Five course corrections are used on 
the subsonic San Francisco to New York City flights and seven on the subsonic 
New York City to London flights. 


Error Models 

i^odels are mathematical representations of the error processes in 
taking measurements. These models are used in the simulation of the real 
Cnoisy) data, which are the inputs to the estimation system, and also as the 
statistical description of the measurement errors used in the Kalman filter 
state estimation computations. In many cases, accurate error models are dif- 
ficult to obtain because it is usually necessary to select a representative 
device and contact the manufacturer for information. Often the manufacturer 
does not have complete information, especially if the device is new, and there 
IS no alternative but to assume a model based only on what little data may be 
available. ^ 


The error models for the various subsystems are summarized in table 2 and 
discussed in the following text along with other models used in the simulation. 

The sensors are sampled and processed by the Kalman filter state estima- 
tor; each sensor considered here provides a measurement of two quantities 
which are sampled simultaneously and called a data set. If the errors in the 
measured quantities have a bias error, it is assumed to be known to sufficient 
accuracy and subtracted from the particular measurement so that, in effect, 
the measurement errors can be considered to be Gaussian distributed random 
variables with a zero mean and a variance of In order to process the 

measurements, a 2 x 2 covariance matrix of the measurement errors, Qit]/) , must 
e given for the particular sensor being sampled. This covariance matrix is 


Qity) 




0120102 o| 


where oi, 02 , and Pi 2 are the standard deviations and correlation, respec- 
ive y, o the errors in the measurements of the two quantities in a data set 


The mepurement subsystem error models are dependent upon assumptions 
regarding their opeptional characteristics. On-board equipment is assumed to 
e modern lightweight digital equipment using modern noise rejection tech- 

forvORTAr/???AN® and accuracy. The ground equipment 

or VORTAC/TACAN stations used in this study are assumed to have a useful 

altitude of 21.336 km (70,000 ft). The Loran-C ground 

assumed to have equipment similar to that in more 
recently established installations. 
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TABLE 2 . — SUMMARY OF SUBSYSTEM ERROR MODELS 


Subsystem 

Error model expression*^ Standard Deviations 

Loran-C 

Latitude error = f (prediction errors + instru- 
mentation errors + coordinate 
converter errors + station 
geometry errors) 

Longitude error = f (prediction errors + instru- 
mentation errors + coordinate 
converter errors + station 
geometry errors) 

VORTAC/TACAN 

Range error = transmitter = 0.55 km 

receiver = 0.055 + 0.165x10 ^ 
Bearing error = 1.3865° x range (km) 

Doppler radar 

Velocity magnitude error = 0.002 |V^[ km 
Drift angle error = 0.2° 

Radio altimeter 

Altitude error = 0.030 km = 100 ft 
Altitude rate error = 0.5 + 0.005|^^| km/hr 

Air data 

Altitude error = 0.183 km = 600 ft 
Altitude rate error = 1.0 + 0.l|h^| km/hr 

Vertical and heading 
reference 

Short period roll, pitch, and heading errors = 
0.25° 

Long period pitch and roll errors = 0.25° or 
0.50° 

Long period heading errors = (0.25° or 0.50°) + 
0.1°/hr 

Body-mounted 

accelerometer 

Acceleration error = 0.01 g 

Altitude and attitude 
control system 

Altitude control error = 0.183 km = 600 ft 
Roll control error = 2.0° 

Heading control error = limited to < ±/2>0° 

Speed control 

Error in speed control = 5,0 km/hr 


Except as noted, all values are 1 o and based on a Gaussian distribution 
with a zero mean. 
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Lovan-C- Loran-C ground transmitting station chains are arranged in geo- 
metric patterns consisting of combinations of one or more three-station triads 
designed to give maximum coverage to a particular area. Each triad of sta- 
tions consists of two slave stations and a common master station and provides 
a position fix in its area of coverage. To conserve equipment, many complex 
station arrangements may be made in which a single station provides multiple 
functions, such as a double slave CDS), master-slave CMS), and so forth, in 
order to provide a large number of position-fixing triads. 


Each master-slave pair produces synchronized pulses. At the receiver, 
the time difference between the arrival of two pulses Csay TD^) from one 
master- slave pair places the receiver on a line of constant time difference 
Ch^erbola) called a line-of-position {LOP) . Similarly, the time difference 
between the arrival of pulses from the same master and a second slave Csay TAj) 
places the receiver on a second LOP. The intersection of these two LOP's is^ 
the position of the receiver. The latitude and longitude of the receiver is 
computed by a coordinate converter from two such time differences measured 
simultaneously and from position coordinates of the master and its slave sta- 
tions. Maximum accuracy is obtained when the LOP' s crossing angle is 90* and 
when the groundwave propagation mode is used rather than the skywave mode, for 
which the signals are reflected from the ionosphere. 

The accuracy of Loran-C position fixes using the groundwave propagation 
mode are dependent upon the following mutually independent error processes: 

Ca) prediction errors in the groundwave propagation; Cb) instrumentation 
errors in the ground station and in the on-board receiver; Cc) coordinate con- 
verter errors; and (d) errors due to the geometrical arrangement of the master 
^nd slave stations and the aircraft* 

Prediction errors arise from assumptions which must be made regarding 

conductivity and dielectric constants of the various media 
oyer which the groundwave must travel, uncertain earth shape, and aircraft 
altitude. Most of these errors can be reduced by calibration and/or compen- 
tion in the on-board computer. Therefore, in this report, the following 
uncertainties (1 o) will be assumed for the prediction error sources: 

1 . Altitude - negligible after compensation based on estimated 
altitude; and 

2. Earth shape and calibration uncertainty overwater and overland - 
0.1 usee (ref. 2 gives a value of 0.06 psec over accurately 
surveyed land) . 

errors at the ground transmitting stations are due to 
0 master and slave stations. Reference 2 gives 

^ this value. The on-board receiver is assumed to be able 

identified''bf^hfn,T to select a particular station triad 

stS^ini ^ tL repetition rate of the master station and two slave 

srar?o*nrISn.r''®^y-"‘ ‘"'"''tinuously tracks the master and two slave stations 
f? as to produce continuously two time differences. The uncertainty in the 
ime di ference measurements of the receiver depends on the signal-to-noise 
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ratio (SNR) of the received signals, which in turn depends on the distance 
from the transmitting station and the ambient noise factor. Following refer- 
ence 2, the SNR is taken as 


10 - 


SNR. = 




50 {d - 130oVb 

max . / 


greater than 10 dB 


1300 < d < 2800 km 
max 


d < 1300 km 
max 


i = (8) 


Here, x and y distinguish the two slave stations of the Loran-C triad and 
^maxv is the larger of the distances (in km) to slave (i) or the master sta- 
tion. The SNR decreases linearly for distances greater than 1300 km and 
reaches the limit of receiver sensitivity (-20 dB) at 2800 km. This estimate 
of SNR is conservative and corresponds to an ambient noise factor of 40 dB 
above 1 yV, a 300 kW transmitter, and propagation across fair soil. 

Following reference 3, the rms receiver errors in measuring TD^ and TD 
can be given by ^ 

where (oi)g '^7 is found from the expressions in the table below, in which SNR^ 
is the minimum master, slave (i) signal-to-noise ratio. 


SNR.; i. = x,y 


i - x^y 

> 10 dB 

0.05 



+0.1957x10-6 

iSNR^)'^ 


-0.8338x10'6 

{SNR^-) 3 

< 10' dB 

+0.4069x10"6 

(5M^)2 


-0.6128x10"2 

(5i?i?^) 


+0.7730x10-1 



The variance of errors in the measurement of TDy; and TDy at the receiver 
due to prediction, ground station, and receiver error sources are then: 

a/ = ( 0 . 1)2 ^ ( 0 . 02)2 + (a^)2^^ i = x,y (10) 

The coordinate converter computes latitude and longitude from the time 
differences TDx ^rid TDy. The conversion is subject to some computation errors. 
Reference 4 gives some examples of conversion errors; if we assume they are 
typical, then the standard deviation of the latitude and longitude conversion 
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errors that could be expected are about 0.46" of arc C2.23xl0'^ rad) and 
0.59" of arc (2. 86x10“® rad), respectively. 


The geometry of the station triad and the aircraft affects the uncertain- 
ties with which the two time differences TD^ and TVy are converted to latitude 
and longitude. The geometrical effect is treated in appendix A where the 
covariance matrix, is derived; this matrix relates latitude and longitude 

uncertainties to the total time difference variances, and and errors 

due to station-aircraft geometry. A final covariance matrix, , that is 

used in the Kalman filtering of the Loran-C data sets can now be given: 


■ %>. * 


(2.23x10“®)2 0 


0 C 2 . 86 x 10'®)2 


( 11 ) 


The maximum usable range of Loran-C is assumed to be 2800 km. At this 
range, the signal-to-noise ratio is at the limit of the receiver sensitivity 
(-20 dB) . In the event that no station triad has a maximum range less than 
2800 bn, the computer skips the position fix. In addition, if the line-of- 
position CLOP) crossing angle, as defined in appendix A, is less than 25® the 
computer attempts to use a more distant station triad. 


VORTAC/TACAN- The military system TACAN operates in the UHF band. From 
the transmitted signals, an aircraft receiver is able to measure its magnetic 
bearing and slant range to the selected station. The bearing is measured by 
phase comparisons using coarse and fine frequencies of 15 and 135 Hz, respec- 
.c slant range is obtained from the Distance Measuring Equipment 
(DME) of TACAN, which measures the transit time of pulses sent from the air- 
craft to a ground transponder and return. 


VORTAC station is the colocation of a civil VHF Omnidirectional Range 
(VOR) station operating in the VHF band and a TACAN station operating in the 
UHF band. The VOR provides the same bearing information as does TACAN, but 
instead measures the phase difference between a 30 Hz signal whose phase is 
varied with magnetic heading and a reference 30 Hz signal. Slant range for 
the civil user is provided by the DME portion of TACAN. 


Twenty-two VORTAC/TACAN stations which lie along the great-circle route 
from San Francisco to New York are used. The physical location of each of 
these stations is taken from a DOD Flight Information Publication (ref. 5) and 
each station is assumed to have either VORTAC or TACAN equipment. Station 
locations are listed in table 3. For this study it is assumed that all sta- 

level and have sufficient power output to be useful 
at 21.336 km (70,000 ft) altitude at a range of 275 km. 


The on-board VORTAC/TACAN receiver is assumed to be a modem digital unit 
capa e o rapid search, lock-on, and rapid channel switching as instructed by 
the on-board navigation computer. The station selection logic used in this 
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TABLE 3.- VORTAC/TACAN STATION LOCATIONS 


Station 

Latitude, north 

Longitude, west 

San Francisco 

37“38’11" 

122°21'36" 

Stockton 

37“50’01" 

121“10'13" 

Sacramento 

38‘’26*37" 

121®33’02" 

Fallon NAAS 

39‘’24*59" 

118“42'10" 

Wilson Creek 

28“15*01" 

114®23’36" 

Bonneville 

40°43’34" 

113®45'24" 

Provo 

40® 16 '30" 

111®56’23" 

Myton 

40®08’42" 

110®07’37" 

Meeker 

40®04'03" 

107®55'27" 

Denver 

39®51'39" 

104®45’08" 

Sidney 

41®05'04" 

103®06’09" 

Grand Island 

40®59'02" 

98°18'52" 

Omaha 

41®10'02" 

95®44'11" 

Des Moines 

41®26'15" 

93®38'54" 

Moline 

41°19'16" 

90®38'17" 

Joliet 

41®32'42" 

88®19'06" 

Southbend 

41®46'07" 

86®19'06" 

Cleveland 

41®21'29" 

82®09'44" 

Philipsburg 

40°54’58" 

77®59'35" 

Selingsgrove 

40“47'27" 

76®53'04" 

Yardly 

40®15'12" 

74®54’29" 

Kennedy 

40®37'50" 

73®46'22" 


report chooses the closest station. This process always resulted in selecting 
a station much closer than 275 km, yet sufficiently distant so that the air- 
craft was never in the "cone of confusion" above the selected station. 

The same error models were used for both the VORTAC and the TACAN range 
and bearing measurements. These error models give the standard deviation of 
the error in the range and bearing measurements about a zero mean. The zero 
mean results from the assumption that constant receiver bias errors have been 
determined, that ground station bias errors versus bearing are known, and that 
a suitable method of correcting for them has been programmed into the on-board 
computer. 

The rms range and bearing uncertainties due to transmitting station and 
receiver errors were taken from references 6 and 7 and are summarized below. 


Measurement 

Transmitter 

Receiver 

Range (km) 
Bearing 

0.055 km (180 ft) 
1.25® 

0.055 km + 0.16x10"^ X range 
0.6® 


The transmitter and receiver errors are assumed to be independent random 
processes so that the matrix for VORTAC/TACAN range and bearing measure- 
ments taken simultaneously is 
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CO. 055) 2 + CO. 055 + 0.165x10-3 x range) 2 | o 


«cv = 


C12) 


CO. 0242) 2 


The uncertainties in range and bearing are assumed to be uncorrelated 
with each other and uncorrelated in time with their previous values, since the 
bias errors versus bearing and range are assumed to be removed by prior cali- 
bration o£ the ground and on-board equipment and subsequent storage of the 
calibration data in the on-board computer. 


Doppler veloG-Cty and drift angle— The Doppler radar system assumed for 
this report is an on-board radar with four downward -looking pencil beams 
arranged much like the four legs of a draftsman’s stool. The antenna radi- 
ating the four beams is space stabilized to the local vertical so that the 
radiation pattern on the surface below is not disturbed by rotational motions 
of the aircraft. By comparing the Doppler frequency shifts from each of the 
four beams with each other, the magnitude of the ground track velocity and the 
Doppler drift angle are determined continuously. The drift angle is the angle 
between the aircraft heading and the ground track velocity vector. 


Tlie Doppler radar system is assumed capable of operating to an altitude 
of 21.336 km (70,000 ft) with continuous output samples as needed by the digi- 
tal computer. Velocity errors are assumed to be uncorrelated in time when the 
sampling interval is 30 sec or greater. The manufacturer's literature indi- 
cates that present Doppler equipment measures groundspeed with 0.2 percent rms 
uncertainty and drift angle with 0.2* rwjs uncertainty. 


Because of the system mechanization, groundspeed and drift angle measure- 
ment errors are essentially uncorrelated. The Doppler drift angle measurement 
is combined with the current measurement of the aircraft heading to provide a 
measurement of ground track heading, GTH. The heading reference is in error 
by an amount independent of the Doppler drift angle so that the combined error 
in measuring the ground track heading is 

GTH (1 o) = [(0.2)2 + 1/2 deg (13) 

The heading reference error variance, (from table 2) dominates the ground 

track heading standard deviation in equation (13), with the result that this 
measurement contributed little to the knowledge of ground track heading or, 
alternatively, the cross-track velocity. It does, however, provide a useful 
check on the performance of the automatic control system. 


The matrix S(*j^) is computed as 




(0.0021^^1)2 

0 


0 




(0.00349)2 + E(^^~) 


where the units are (km/hr) 2 and (rad) 2, 


(14) 
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Radio altimeter^ A radio altimeter is a downward- looking radar ranging 
system which measures the range from the aircraft to the surface below. The 
recent design trend is to integrate some of the functions of the radio altim- 
eter with those of the Doppler radar. This integration allows a spce stab- 
ilized antenna for the radio altimeter at low cost. The radio altimeter 
assumed for this report is a high power device capable of measuring altitude 
above the local surface, up to 21.336 km (70,000) ft), A radio altimeter was 
assumed only for flights overwater, since terrain effects would not have to be 
considered. The radio altimeter operates continuously and its output is 
sampled under control of the on-board digital computer. Reference 8 gives a 
1 o uncertainty of about 0,0061 km (20 ft) in the measurement of altitude for 
a system which has been adjusted immediately prior to the measurement. Since 
the system assumed here would not normally be adjusted prior to a measurement, 
the vms altitude uncertainty is taken as 0,0305 km (100 ft). 


Altitude rate is obtained from the radio altimeter by measuring the 
Doppler shift in the carrier frequency, which in this application also sug- 
gests combining the functions of the radio altimeter and the Doppler radar. 
The manufacturer’s literature for this type of equipment indicates that the 
uncertainty in the measurement of altitude rate is about 0,5 km/hr plus 
0,5 percent of the altitude rate. Therefore, the matrix used for pro- 

cessing altitude and altitude rate is 


ecV = 


(0.03)2 

0 


0 

(0.5 + 0.005|^j,|)2 


(15) 


The errors in altitude and altitude rate are treated as uncorrelated because 
this relationship is implied in the literature. For example, in some of the 
more modern systems, the modulation of the carrier is sinusoidal and the alti- 
tude is found from the relative phase of the Bessel sidebands and the corre- 
sponding transmitted signal. The altitude rate, on the other hand, is mea- 
sured by the Doppler shift of the carrier frequency. 


Aiv data- An air-data computer which has aerodynamic and thermodynamic 
sensors is used overland to compute such quantities as altitude and altitude 
rate, airspeed, true airspeed, angle of attack, sideslip, and Mach number. 

The discussion at the end of the series of articles in reference 6 indicates 
that with very accurate calibration equipment, altitude from an air-data com- 
puter can be expected to have an vms uncertainty of about 0.183 km (600 ft) 
at 21.336 km at (70,000 ft) altitude. Very little information could be found 
regarding the accuracy of altitude rate from air-data computers. It was arbi- 
trarily assumed that with the various compensations applied by the air-data 
computer, the rms uncertainty in altitude rate is the same at 10.668 km 
(35,000 ft) and 21.336 km (70,000 ft) altitude. The value chosen was l.Okra/hr 
plus 10 percent of the estimated altitude rate. Therefore, the matrix 

for simultaneous processing of altitude and altitude rate was 




(0.183)2 

0 


0 

(1.0 + o.ilhj,|)2^ 


(16) 
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The altitude and altitude rate errors were assumed uncorrelated for simplicity. 
In addition, the uncertainty chosen for the altitude rate is sufficiently 
large so that the system essentially ignores the measurement. The poor mea- 
surement of altitude rate causes no difficulty because the guidance system 
attempts to maintain constant altitude while en route and therefore relies 
mostly on altitude measurements. 

Vertioal and heading referenoes- The aircraft is assumed to be equipped 
with a vertical (or attitude) reference consisting of a vertical gyroscope 
with gravity sensing of the vertical. Vertical reference of this type, using 
third-order leveling and corrections for Coriolis acceleration, are capable of 
relatively good accuracy in pitch and roll. In this report, the vertical ref- 
erence is capable of measuring short period changes (a few seconds) in pitch 
and roll to an ms uncertainty of 0.25*. Over longer periods of time, the ms 
uncertainty is taken as either 0.25* or 0.5*, depending on the quality of the 
vertical reference. The mis errors are assumed to have a Gaussian distribu- 
tion with a zero mean. 


The heading reference can be mechanized in combination with the vertical 
reference, or, for example, with a directional gyroscope with magnetic sensor 
slaving, as has been assumed here. As with the vertical reference, the ms 
uncertainty over short periods of time was taken to be 0.25°; over longer 
periods of time, either 0.25* or 0.5* was assumed, depending on instrument 
quality. Since the magnetic variation is known to an uncertainty of about 
0.1^ over most of the earth's surface, a pseudodrift of the gyroscope of 
0.1 /hr was assumed between course corrections to simulate the effect of this 
uncertainty. 

Body-mounted aaoelerometer- An accelerometer is mounted at the center of 
gravity with its sensitive direction along the aircraft longitudinal axis. 

This accelerometer is used by the automatic control system to control the air- 
craft speed during cruise and to measure speed changes during course correc- 
tions. The accelerometer output is assumed to have a Gaussian-distributed 
error with a zero mean and a standard deviation of 0.01 g. Gravity-induced 
acceleration due to aircraft pitch attitude is compensated for in the on-board 
computer. 

Automatic control system— An automatic control system (autopilot) is 
assumed to be in operation throughout the flights. Its functions are the fol- 
lowing: (1) to maintain the attitude of the aircraft, that is, the pitch, 

roll, and commanded heading; (2) to hold the prescribed altitude; (3) to per- 
rora the required maneuvers to execute the commanded course corrections; and 
(4) to maintain speed or respond to commanded changes in the aircraft speed 
UTing course corrections (autothrottle) . This automatic control system was 
^ permit constant rate climb-out and descent maneuvers even 

t oug t ey were not performed in this study. The course correction computa- 
tions include a pitch change computation for this purpose. 

To simplify computations, the control system on all flights was designed 
o. (1) maintain the cruise altitude to within an error whose standard devi- 
ation is 0.183 km (600 ft); (2) maintain the roll angle at the commanded value 
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to within an error whose distribution is Gaussian with a zero mean and a stan- 
dard deviation o f 2. 0°; (3) maintain the heading such that the heading error 
is limited at ±/2.0° and has a Gaussian distribution with a zero mean and a 
standard deviation of ±/2 . 0* between these limits; and (4) control the speed 
to within an error which is Gaussian-distributed with a zero mean and a stan- 
dard deviation of 5 km/hr. 

Aircraft velocity vector wander between measurement data-sets— During the 
interval between processing two successive measurement data-sets at, say, t\_^ 
and the aircraft velocity vector wanders from that predicted by the equi- 
librium transition matrix derived in appendix B. This wander is due to random 
errors in controlling the effects in velocity of random disturbances in wind, 
air density, temperature, and so forth, which occur during 

The distribution of the wander error from the commanded velocity can, 
therefore be calculated from the automatic control system error model, and the 
acutal wander error is simulated by sampling this distribution. 

Appendix C treats these error sources and derives expressions to modify 
the actual and estimated aircraft states and the E and P covariance matrices 
to account for the effect of wander from the equilibrium state. Once the com- 
putations in appendix C are complete, the Kalman filter computations described 
in the section on State Estimation are started. 

The control system errors are modeled as correlated in time (appendix . 
The numerical value of the correlation decay coefficient, a, was unspecified 
in appendix C. Early in this study it was decided to determine a single 
choice of a for all of the data runs since there was a large number of other 
variables to be considered. This choice was made by evaluating the effect of 
parameter a on system performance. Several runs were made with various values 
of a and it was found that performance was insensitive to values between 30/hr 
and 120/hr. As a result, the value a = 60/hr was chosen as being representative. 

State Estimation 

The objective of the state estimation system is to provide an accurate 
estimate of the actual aircraft state. To accomplish this objective, a Kalman 
filter is employed to process simulated measurement data-sets which are 
derived from the simulated actual aircraft state and statistical descriptions 
of the measurement noise. A block diagram of the computation flow involved in 
the state estimation is shown in figure 4. 

The block diagram is divided into three portions. The upper portion 
illustrates the computation flow for the simulation of the aircraft state 
(actual state), including the computation of the actual add measured data-sets. 
The lower portion of the block diagram illustrates the computational flow in 
the reference path simulator. The center portion illustrates the computational 
flow of the estimation system. The estimation system is composed of an esti- 
mated path simulator (analogous to actual and reference path simulators), the 
Kalman filter, and the computation of the estimated data-set and residual 
error-set, which is used to compute corrections to the estimated state. 
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Figure 4.— Block diagram of ccsnputational flow for state estimation. 

In the following detailed discussion of the computations embodied in the 
simulated sptem, the discussion is centered about the Kalman filter, since 
the simulation of both the actual and reference states involve only simple 
updating when a data-set is processed. The discussion is divided into five 
sections. The first four sections deal with the major divisions of the Kalman 
filter computations and the fifth section deals with computations for esti- 
mated performance evaluation. These sections arei 


1. Updating the covariance matrices and the states from ^ to 

2. Computation of the weighting matrix, ^ 

3. Updating the covariance matrix, P(t7)^ 

4. Computation of the corrected estimated state. 

5. Computation of r and V for performance evaluation. 


(^dating the oovcxriccnoe matrvces and the states from tv^i to ti— The 
covariance matrix of the errors in the prior estimated state, P(tT, the 

SItes^"^rt?^^y^^aSd Iht deviations between the actual and ri^rence 

states, and the three states Xj^, and X^ are available after 
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processing the (fe - 1^ the data-set. It is now desired to process the current 
data-set at time t}>_. To initiate the processing the covariance matrices P and 
Pj the estimated state, the actual state, and the reference state (6 x 1 vectors) 
must be updated through prediction from to the current time This 

process is accomplished by use of their respective transition matrices 

, where t = /i, P, and E. The transition matrices are found from 
equations derived in appendix B. 

The covariance matrix, ^)» is updated in time, using the equation 

P(*^) - 

where the matrix AP (derived in appendix C) accounts for the errors in measur- 
ing the wander of the velocity vector during the time period 

The matrix Rit-0 used for performance evaluation is the covariance matrix 
of the error between the actual and reference states. It is updated using the 
equation 

where AP is a covariance matrix that accounts for the effect of random 
perturbing forces on the actual path of the aircraft between data-set samples. 

The actual, estimated, and reference states are updated by their respec- 
tive transition matrices using the following equations: 




(19) 



(20) 



(21) 


where equations (19) and (20) also indicate the addition of the vectors hX/[ 
and AZgp to produce the updated actual and estimated states. The vectors AZ^ 
and are the actual and measured deviations from the paths predicted by the 
transition matrices ^E^^ky'^k-O > respectively, during the time. 

They are used to account for the random perturbing forces on the 
aircraft which are not modeled by the transition matrices and have covariance 
matrices AP and AP, respectively. 

Computation of the weighting matrix, K(ti^)— Mt&r the updating computa- 
tions are complete, the Kalman weighting matrix, K{t]f) , is computed as follows. 
First, the measurement Jacobian matrix, H{t-0 , is evaluated from partial 
derivatives of each of the measurement quantities in the particular measure- 
ment data-set being sampled. The P matrix has two rows since each measurement 
data-set consists of two measurement quantities. Expressions for the R{t}f) 
matrices are derived in appendix D for each (Loran-C, VORTAC/TACAN, etc.) mea- 
surement data-set. Next the elements of a covariance matrix of the measurement 


21 



errors, Qit]^ , (derived in the Error Models section) is evaluated for the 
appropriate data-set. Finally, the Kalman weighting matrix, K{t.) , is com- 
puted from the expression: ^ 

ifCV = P(tp/(y (22) 

where the indicated matrix inversion involves only a 2 x 2 matrix. 

Updating the oovariance matrixy F(tj^)- After K{t-0 has been computed, 
P(*fe) is computed to reflect the reduction of the uncertainty between the 
actual and estimated states as a result of the information gained from the 
measurement data-set. The computation is carried out as follows: 

P(tj^) = [J - ( 23 ) 

where is the corrected matrix. After a delay to the time when the next 

data-set sample is to be processed, P{tp) will become P{t]^.i) , and is used in 
equation (17) as the initial information for the same sequence of updating and 
processing computations applied to the next measured data-set. 

Computation of the oovreoted estimated state— The estimated state is 
updated to reflect information gained from processing a measured data-set. To 
simulate a measured data-set, an actual data-set is computed from the geometry 
of the aircraft relative to the earth and the actual state variables. Then 
measurement noise, which is generated according to the variance given in Q{tif\ , 
is added to the actual data-set to produce the measured data-set, y. This 
process differs from the procedure which would be followed on-board, since the 
actual state is not known and the measured data-set would be obtained from 
on-board equipment. The estimated data-set, y, is computed from the same geom- 
etry of the aircraft with respect to the earth as was used above and the esti- 
mated state variables. A residual error-set, y - y, is then computed by 
taking the difference between measured and estimated data-sets. The estimated 
state is updated using the following computation: 

- y) (24) 

where is the updated estimated state. 

The text above has given the detailed steps involved in processing a data- 
set with the Kalman filter to produce an updated estimate of the aircraft 
state. The difference between the actual state and the estimated state, often 
called the estimation error, may be considered a measure of estimation system 
performance since a perfect estimate of the aircraft state would result in 
zero estimation error. 

Computation of r and v fov peipfoTmanae evaVuxxtion— Figure 4 indicates the 
computation of the estimation error vector, X{t-0 , (6x1 vector) which is 
obtained by subtracting the estimated state from the actual state. The quan- 
tities r and v (3 x i vectors) are computed by taking the square root of^the 
sum of the squares of the position and velocity terms, respectively, of X{tf) . 

K 
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This description completes the computational flow illustrated in figure 4 
and is the sequence required to process a simulated measurement data-set; this 
process is repeated for each new measurement data-set. 


Flight Path Control 

After repeated use of the estimation procedure described above, a reason- 
ably good estimate of the aircraft state should be available from the digital 
computer; the aircraft speed and heading can then be corrected using the esti- 
mated state and the known reference state. In this report, it is desired that 
corrections to the flight path return the aircraft to the reference path and 
maintain the reference position and velocity in order to arrive at the desti- 
nation at a fixed time. 


One way to accomplish this goal is by using the two-maneuver correction 
scheme shown in figure 3. When the aircraft arrives at the scheduled time of 
the first maneuver, tj, the fixed time increment, ht\, and the estimated state 
are used to compute a transition matrix, from to as illus- 

trated in figure 5. The fixed-time-of-arrival correction equations (see appen- 
dix E) are then used to compute Av from the estimated and reference states. 


At I At 2 
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Figure 5,— Block diagram of fixed-time-of-arrival flight path control. 

The increment Au is the commanded change in the velocity vector required for 
the aircraft to arrive at the desired intercept point on the reference trajec- 
tory at time ta. 














This commanded velocity correction is transformed into pitch, heading, 
and speed corrections and compared to the execution uncertainty for each of 
these three variables. The execution uncertainty is considered to be the min- 
imum response the control system can make. In pitch and heading, it is the 
short-term rms control error (0.25*) and in speed, it is the rms speed control 
error (5 km/hr) . If any one of the commanded changes is less than its respec- 
tive execution uncertainty, that commanded change is cancelled and only the 
remaining commands are executed according to the course correction implemen- 
tation model (appendix E) . The resulting changes are added to the actual and 
estimated states to implement the course correction. 

After the course correction, F(t?c) and F(-tfe) become P(t^) and and 

must be adjusted for the effects of a course correction. This method is 
detailed in appendix E, but will be described in general here in order to com- 
plete the description of the block diagram in figure 5. 

Following a course correction, the estimation errors are increased as a 
result of errors in measuring the executed course correction. Accordingly, 
P(tp is changed, using: 




fo 0 -] 



L 


where P(t;;.) is the new covariance matrix. The 3 x 3 matrix P(nn^) is the 
covariance matrix (appendix E) of errors in measuring the executed correction. 
Next, to update the P(t^), a correction matrix, G, foniied when Ay was computed 
in the fixed-time-of- arrival equations, is used along with P(t, ) and as 

follows: ^ ^ 

= (J + (7)[P(t“) - P(tp](J + + P(t^) (26) 

where P(*p is the new matrix which replaces P(i^. 

At time t 2 , the aircraft will be near the desired reference path and the 
second maneuver of the two-maneuver course correction sequence must be imple- 
mented to aline the aircraft along this path. The second maneuver is executed 
in the same manner as the first, using the same flow of computations as shown 
in figure 5, with the following exceptions: 

1. The intercept point is the desired destination. 

2. The time-to-go to the desired destination replaces At, 

3. The transition matrix is $„(t,, t ), where t, is the reference time at 

the destination. a 2 a 
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Initialization o£ the Reference^ Estimated, and Actual States 


The initial reference state is found from the desired cruise altitude, 
the desired reference velocity magnitude, and the initial and final positions; 
it serves as a basis for the determination of the initial estimated and actual 
states. The initial reference position and velocity are derived from the posi-* 
tion of the starting point, the direction of a great-circle route from this 
point to the destination, and the cruise speed of the aircraft. The computa- 
tions are done in the following way: First, unit vectors from the center of 

the earth to the starting point, and the final position (destination), ^ 2 ^ 
are found from their latitude and longitude. Second, the cruise initial posi- 
tion vector, L^, is computed from 

Lr" 

where h is the desired cruise altitude and Eg is the radius of the earth. 

Third, a unit vector, u, perpendicular to the plane of J-i and )l2> computed 
from: 


u 


Jlj X ^2 


1^1 ^ ^2 


(28) 


Fourth, the reference angular velocity vector (with respect to the center of 
the earth) is then 


where |7/?| is the desired magnitude of the reference groundspeed. Finally, 
the initial reference velocity vector is then 

7„ = 0) X L (30) 

ti ri 

The best estimated state at t = 0 is assumed to be the initial reference 
state. As a result, deviations between the actual and reference states, and 
the actual and estimated states, are initially identical. These initial devi- 
ations are assumed to be random variables with zero means and covariance 
matrices i?(0) and P(0), respectively, which are also identical as a result of 
the above assumption. 

The matrix i?(0) if found by assuming that the Doppler radar and on-board 
digital computer have a secondary "navigate" mode in which the along-track and 
cross -track components of velocity and computed position are use by the con- 
trol system to reduce the effect of random disturbances such as wind and so 
forth, and to maintain control over the velocity vector from t = -10 min to 
t = 0, where navigation and guidance is taken over by the primary system. At 
this time, i?(0) is assumed to be the sum of a diagonal matrix representing the 
errors in initializing the secondary navigation mode and a second matrix 
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typical of the AR matrices, computed between measurement data-sets and repre- 
senting the errors in the secondary navigation system from t = -lo min to t=0. 
This computation was done by choosing 6R as a diagonal matrix, which, when par- 
titioned into 3x3 submatrices, is given by 

ri 0 1 

= (31) 

Lo 100 j 

where SR is in the ALT, ALONG-TRACK, CROSS-TRACK coordinate system. The 
second matrix, Ai?(0,-10), was generated by evaluating the AR matrix (appen- 
dix C) for an attitude uncertainty of either 0.50° or 0.25° from t = -10 min 
to t = 0 min. The desired R(_0) matrix is then found from 

i?(0) = Ai?(0,-10) + SR (32) 

The initial actual state deviations from the reference state are gener- 
ated from i?(0) using a random number generator; these samples are added to the 
reference state to form the actual initial state. The initial deviation used 
tor the sample cases of this study are summarized in table 4. 

TABLE 4.- INITIAL DEVIATIONS OF ACTUAL STATE FROM REFERENCE STATE 


r ' " 

Attitude 

uncertainty, 

deg 

Altitude 

Along-track 

Cross-track 

L Supersonic 


t^osition. 

0.50 

1.353 

0.209 

-1.732 

km 

.25 

.722 

.208 

.924 

Velocity, 

.50 

-9.065 

1.859 

29.025 

km/hr 

.25 

-5.281 

1.924 

19.000 


Subsonic 


Position, 

.50 

.556 

.208 

-.712 

km 

.25 

.376 

.208 

-.481 

Velocity 

.50 

^ -4.370 

1.922 

16.592 

km/hr 

.25 

-3.481 

1.923 

14.239 


. be noted that the four cases in table 4 are essentially one 

individual case with identical initial conditions at t = -10 min. The differ- 
ences at t = 0 are due to the differences in velocity and the attitude 
uncertainty . 


PERFORMANCE EVALUATION 


Performance evaluation computations were included as an integral portion 

programmed in FORTRAN IV and executed on an 
result of past experience with state estimation algo- 
p ograjnmed on this computer, some special precautions were taken to 
preven round-off errors from dominating both the simulation results and the 
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performance evaluation computations. These precautions included double pre- 
cision vector dot and cross product calculations, forced symmetry of the 
covariance matrices, and forming the elements of a matrix product by accumu- 
lating products according to the sign prior to the final addition to form the 
elements. In addition, performance calculations were done in double precision. 
All other computations were done in single precision. 

For on-board computation, it seems likely that computers with shorter 
word lengths than the IBM 7094 would require at least partial control of numer- 
ical round-off in the vector and matrix operations. In addition, a method of 
preserving symmetry of the P and R matrices is recommended. Recent advances 
with computational techniques have made the "square root" form of the Kalman 
filter appear attractive for on-board computations. Although some double pre- 
cision computations may be required with this method, the symmetry, for exam- 
ple, is assured. 


Estimated and Actual Flight Paths 

Three flight paths were computed: a reference path, an estimated path, 

and an actual path. The reference flight path is the desired path of the air- 
craft. The reference state variables can be computed for any elapsed time^ 
along this path and are therefore considered to be known variables with which 
the estimated and actual values can be compared. The estimated flight path 
was computed and updated in position and velocity by use of a Kalman filter 
which processed discrete measured data-sets from the various digital subsys-^ 
terns as selected by the on-board computer. After each measurement data-set is 
processed, the resulting new estimated state is the systems* estimate of the 
actual state. New estimates are constantly required since both the actual and 
estimated states are perturbed from the states predicted by their respective 
transition matrices in the period between measurement data-sets and by the 
execution of course corrections. 

Performance of the state estimation and course correction systems are 
evaluated by comparing the estimated state with the actual state and by com- 
paring the actual state with the reference state. Comparisons are made indi- 
vidually and statistically. Individual comparisons between the estimated ^ 
state and the actual state are the errors in the estimate and are denoted r 
and V. The statistical errors in the estimate are denoted '^rmS’ 

Individual comparison between the actual state and the reference state vari- 
ables are denoted r and v and the statistical comparisons are denoted and 


Statistical Information for Performance Evaluation 

The statistical information for performance evaluation of the state esti- 
mation system can be derived from the covariance matrix P. The statistical 
information for performance evaluation of the combined state estimation, guid- 
ance and control systems is provided by the statistics of the actual path 
deviations from the reference path contained in the covariance matrix R. 
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These two matrices are 6^6 and the statistical data is obtained from them by 
first partitioning into 3 x 3 submatrices as follows: 




-Pl P2- 




P = 

.^3 


(33) 

and 

R = 

“Pl R2~ 
L^3 ^4_ 


(34) 

From the P submatrices we 

obtain the following: 




(35) 


^rms " 

(36) 

and from the E submatrices 

we obtain in a similar fashion 



2 * = (trace 

rms ^ 

(37) 


^rms " 

(38) 


These quantities represent the total statistical position and velocity errors. 
In the case o£ and a statistical measure of the total error in the 

estimated position and velocity is obtained. In the case of and Vrme. a 

statistical measure of the total actual deviations from the reference state 
are obtained. At this point, it should be made clear that the statistical Tor 
ensemble average] performance of the system is obtained from a single run on 
the computer and that, should a large number of computer runs be made, it 
would be expected that ensemble averages of the individual deviations r, u 


2^3 and V would be given by i* ^ v 

^ rms 


rms 


VJ 

rms^ 


^rid respectively. 


Course Control Performance 

Evaluation of course control performance is based on the ability of the 
system to control position and time errors continuously along the entire ref- 
For ease in plotting and computation, the lengths of the refer- 
ence paths were adjusted slightly so that the time en route would be a more 
convenient value. As a result, the en route reference times were as follows: 

tm^don.^SOO rtl?’"' aircraft. New Yeri City 


28 



Special scales are provided on the figures for evaluation of en route 
timing performance along the entire reference flight path. These scales allow 
determination of the en route timing error at any point in time along the ref- 
erence path by conversion of t*, rrmsa ^.nd r into equivalent time errors. 

The en route timing errors are denoted as respec- 

tively, and are computed from: 


6* 

kr 

(39) 

rms 

rms 


6* = 

kr 

(40) 

6* 

kr 

(41) 

rms 

rms 

11 

kr 

(42) 


where the constant k for the supersonic aircraft is equal to 1.2 sec/km and 
for the subsonic aircraft, 3.6 sec/km. Each of these time error quantities is 
the time for the aircraft to travel the total distance given by the respective 
position error quantity at the reference velocity. In particular, the quan- 
tity 6*27778 time required to travel from the rms actual position to the 

reference position at a specific time en route. The quantity 6t is the actual 
en route time error for one member of the ensemble of times represented by 
6*27778 • The quantity 6*27778 is the uncertainty in the knowledge of the time 
necessary to travel from the actual position to the reference position and 
means that the en route time error would be known on board the aircraft with 
an uncertainty given by The quantity 6* is one member of the ensemble 

of errors represented by 

DISCUSSION OF RESULTS 


The navigation system performance of subsonic and supersonic aircraft was 
evaluated for four cases. Flights were routed over a great-circle route over- 
land from San Francisco to New York City and a great-circle route overwater 
from New York City to London. Attitude uncertainties (whose major effect is 
on heading control) of 0.50^’ and 0.25® were used in conjunction with a fixed- 
time-of-arrival (FTA) course control system. System performance of the four 
cases will be discussed in the following six areas: (1) state estimation per- 

formance; (2) effect of VORTAC/TACAN data-set samples; (3) effect of course 
corrections on the rms velocity estimation error; (4) individual case estima- 
tion performance; (5) control of the actual state; and (6) en route 
timing errors. 


San Francisco to New York City — Subsonic Cruise 

In this case, the aircraft reference velocity and altitude are 1,000 km/hr 
and 10.668 km (35,000 ft). The data-sets processed by the Kalman filter were 
from VORTAC/TACAN, Doppler radar, and barometric altitude information from an 
air-data computer. These data-sets were taken every 3 min. 


29 



state estimation performance-- Control of the actual aircraft state for an 
effective FTA system is dependent upon the accuracy of estimating the aircraft 
state at the time the course corrections are executed. State estimation per- 
formance with two different attitude uncertainties (0.50® and 0.25®) is shown 
in figure 6. The figure shows an individual total position estimation error 


Subsonic aircraft 

San Francisco to New York City 

0.50 deg 

0,25 deg 


Type of data 
sampled 


Vortac 

/^Doppler 


« • • • 


"Altitude 


HccK Hcch Heel- Hech Hech 



Figure 6. Comparisons of r and for A.U. = 0.25 and 0.50°. 
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for the individual ensemble member> i* (top), and the statistical total urns 
position and velocity estimation errors, ^rms position 

estimation error for an attitude uncertainty of 0,50 varies from about 0.60 
km to about 1.7 km; for an attitude uncertainty of 0.25”, it varies from 0.5 
km to 1.2 km. 

Total rms velocity estimation errors (middle) are nominally of the order 
of 10 km/hr, with variations from 9 km/hr to 15 km/hr for an attitude uncer- 
tainty of 0.50”; for an attitude uncertainty of 0.25”, the total rms velocity 
estimation errors are nominally about 6.5 km/hr with variations from about 
5 km/hr to about 10 km/hr. 

Effect of VORTAC/TACAN data-set samples- The data presented in all the 
figures are what existed just after a data-set was processed by the Kalman 
filter and a new estimate of the aircraft state was obtained. Dots at the 
tops of certain figures are used to indicate the type and timing of data-set 
processing. For example, a dot at the coordinates of 30 min time en route 
and in the VOR row would mean that a VORTAC/TACAN data-set was processed at 
30 min time en route. This allows a determination from the figure of the 
effect on r , y , and r of each type of data-set at any time en route 
along the r^irence^path . 

The VORTAC/TACAN data-set samples are the primary source of position in- 
formation in the ALONG-TRACK and CROSS-TRACK directions. During the initial 
settling period of the estimation systej^, the large reductions in r and Vrms 
occur at VORTAC/TACAN samples. Later, is held in the range of 0.5 to 

1.5 km largely due to the reduction of r^^g to minimum levels at each VORTAC/ 
TACAN sample. For the individual case, ^ can differ from this trend due to 
unfortunate combinations of random errors, but as seen in figure 6, most of 
the minimum points of r occur at VORTAC/TACAN samples. 

All of these curves exhibit the ragged appearance characteristic of dis- 
crete estimation systems. Amplitude details can be varied by the frequency 
of processing and by the data type used. These curves do not exhibit a step 
reduction in estimation error at the time each data-set is processed because 
only the estimation error existing after the data-set was processed is shown. 
This method of presenting results has the disadvantage of not showing explic- 
itly the improvement made by the processing of each data-set, but experience 
has shown that the performance after processing a data-set is of far greater 
interest and the clarity of the figures is greatly enhanced. 

Effect of course corrections on the rms velocity estimation error— In 
figure 6, the increase in the rms velocity estimation error, Vrms> due to 
uncertainties in the execution of a course correction is shown by a peak in 
the curve just after the second maneuver of each course correction. This 
peak is most visible on the curve for an attitude uncertainty of 0.25”. Peaks 
at the same time can be seen, though less clearly, in the curve for 0.50”. 
Peaks also occur after the first maneuver of each course correction. But in 
this case, a Doppler radar data-set was processed by the Kalman filter imme- 
diately following the execution of each first maneuver and effectively removed 
the velocity uncertainty introduced by the course correction; thus the 
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velocity uncertainty was reduced to approximately the same level existing 
before the execution of the first maneuver. 

Individual case estimation^erformance— The total position estimation 
error for the individual case, r», is the vector sum of the three error com- 
ponents of the individual position estimation. The individual estimation 
errors influenced the individual course corrections which were executed on 
this particular flight. In particular, they are estimation errors which ex- 
isted at the time of the computation of the individual course correction 
maneuvers and relate directly to the system's performance in controlling the 
actual state and the corresponding timing errors. 

The individual velocity estimation errors are also involved in the course 
correction computations but, since their influence is relatively minor com- 
pared to that of the position estimation errors, the curves are not shown. 

Control of the actual state— The rms position estimation error, '^xms> for 
an attitude uncertainty of 0.50° is compared in figure 7 (bottom curve) to the 


O Total position error (Ti-rris) fi'om the great circle reference course, 
standard deviation (I<t) 

□ Along track component of Trms 

^ Cross-track component of Trms 

O Total error In the estimation of position (7rms) with respect 
to the actual position 


Data sampled: Vortac/Tacan,doppler radar^ barometric altitude 
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Figure 7 rms errors in position and position estimation subsonic 
case (1000 km/hr). San Francisco to New York City A.U. = 0.5°. 


actual position deviation from the reference position, (top curve) and 

the along-track and cross-track components of Since the vertical com- 
ponent of small, the total is very nearly the vector sum of the along- 

track and cross- track components. Data points in figure 7 are plotted at^ 
about 6 to 9 min intervals and show that even though the position estimation 
error is only about 1 km, the minimum actual error, is about 4 km. 

effect occurs at the intercept point (see fig. 3), which occurs at the time o 
the second maneuver of each course correction and shows that with five course 
corrections, position estimation performance is about four times better than 
position control of the actual state. It might be expected that r^rns would 
more closely approach this time, but it does not for three reasons. 

(1) the course correction computations are based on estimated state informa- 
tion; (2) course corrections are not executed perfectly; and (3) errors in the 
control system will allow external forces to deviate the aircraft from the 
computed path to the intercept point. These three error sources cause posi- 
tion errors which, to a first order approximation, increase linearly with 
time. This effect suggests that position control would be improved by in- 
creasing the frequency of course corrections and that the optimum performance 
(statistically) could be achieved by an adaptive system, in other words, by 
making a course correction following the processing of each data-set, pro- 
viding the course correction threshold requirements were met. 

It is of theoretical interest to pursue this line of reasoning further. 
Suppose there were no errors in the course correction computations due to un- 
certainties in the state estimates; then for conditions cor;£esponding to those 
of figure 7, it would be expected (statistically) that for ^ 

rms position error would have a lower limit given by 

r = r + |7o|At(0.5/57.3) 

= 1 + (1000) (6/60) (0.5/57.3) = 1.9 km 

This position error is based upon the premise that a two-maneuver course cor- 
rection is made every 6 min. As a practical matter, experience has shown that 
with the present course correction threshold, course corrections would be made 
much less frequently than every 6 min. As a result, the 1.9 km value would 
not be approached. Even so, a value substantially less than 4 km is feasible 
through a higher frequency of course corrections, such as would be provided 
through an adaptive system. The effects of increased frequency of course 
corrections was not investigated explicitly in this study, since it is clear 
from the foregoing discussion that a significant improvement in performance 
can be obtained by this means; investigation of this subject was left for 
further study. Instead, the use of a fixed near-minimum number of course cor- 
rections shows much more dramatically the individual effects of the various 
components of the system and the data-set sampling schedules. 

The two horizontal lines in figure 7 at 7.4 km (4 n. mi.) and 18.5 km 
(10 n. mi.) are present separation standards (cross-track) within the United 
States (ref. 9). When the aircraft has distance measuring equipment, the 
7.4 km (4 n. mi.) standard applies. The cross- track component of meets 

this separation standard on a 1-a basis. 
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The quantity r>^^, which is indicative of actual state control, is com- 
pared for attitude uncertainties of 0.25" and 0.50" in figure 8 (bottom curve). 

Subsonic aircraft 

San Francisco to New York City 

0.50 deg 

0.25 deg 

, . Vortac 



Doppler 

Altitude 

57.6 

54.0 

50.4 

46.8 

43.2 u 

tfj 

39.6 

o 

36.0 t: 
o> 

32.4 

c 

28.8 E 

25.2 Z 

21.6 o 

18.0 " 

14.4 « 

10.8 

7.2 
3.6 



Figure 8.- Comparison of and r for A.U. = 0.25 and 0.50“. 

Also shown is the total individual case actual position deviation from the 
reference, r, for each of the two attitude uncertainties. Figures such as 
this one are useful in judging the effectiveness of a particular number of 

meeting specific cross-track error requirements. In 

actual nairfrL''tr’ deviation of the 

for ^ reference path varied from about 3.8 km to about 8.3 km 

tor 0.50 altitude uncertainty and from about 2.8 km to about 6.8 km for the 


34 


1 



0.25® case. The upper curve, v, (individual case) is of interest because it 
dramatically demonstrates the wide variation in actual flight path deviations 
with a given control system accuracy and a small number of course corrections. 
It also calls attention to the discrepancy in the statistical (rms) perform- 
ance and the actual performance which can occur; the rms plot indicates quite 
adequate control, whereas in the individual case, control errors were signi- 
ficantly larger. This effect is due primarily to unfortunate combinations of 
random errors occurring in the control system. Better control of the indi- 
vidual case could be obtained by: (a) using an adaptive correction schedule 

which would make a course correction whenever the groundspeed or heading 
change required exceeded a given threshold; or (b) using the existing fixed 
schedule and imposing an upper threshold on the estimated deviation which, 
when exceeded, would cause an additional correction to be implemented. 

En route timing errors^ The course correction system used in this report 
is intended to reduce en route timing errors at all points along the refer- 
ence path to acceptable values. Figures 6, 7, and 8 have scales on them so 
that en route timing errors may be found at any point along the reference 
path. The definition of en route timing error used herein is the time it 
takes to travel from the actual aircraft position to the reference position 
at a particular time and at the reference velocity. 

The approximate range of the rms en route timing errors, 

OT8 uncertainty in these errors, are given in the table below: 


Attitude uncertainty. 

St 3 sec 

6^ 

deg 

■ms’ 

rms sec 

0.50 

13.5 - 30.0 

1.8 - 6.0 

0.25 

10.5 - 24.5 

1.8 - 4.1 


The table shows that the major effect of the attitude uncertainty is on 
and results from decreasing control system performance with increasing atti- 
tude uncertainty. 


San Francisco to New York City - Supersonic Cruise 

In this case, the aircraft is supersonic with a reference velocity of 
3,000 km/hr. The flight path is from San Francisco to New York City along a 
great-circle path at an altitude of 21.536 km (70,000 ft). 

State estimation performance— In figure 9, r^^, ^ shown 

for attitude uncertainties of 0.50® and 0.25®. The data-set sample rate is 
1/min, 
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rms position estimation 





and 


The approximate ranges of the rtns position estimation error, 
the velocity estimation error, shown in the following table. 


Attitude uncertainty, 
deg 

r , km 

rms 

y , km/sec 

rms 

0.50 

0.35 - 1.25 

19.5 - 26.5 

0.25 

0.35 - 0.90 

9.0 - 18.5 


As can be seen from this table, the effect of the attitude uncertainty is pri 
marily on 

Effect of VORTAC/TACAN dataset samples on state estimation- When figure 
9 is compared to figure 6, which shows the corresponding data for the subsonic 
case, it is found that in general, the effects of VORTAC/TACAN data-set sam-^ 
pies on '^rms ^-re much the same as in figure 6. In addition, there is 

a reduction of over that of the subsonic case. This reduction is due to 

the shorter time intervals between data-set samples, which lessen the effects 
of external disturbances since velocity and acceleration disturbances have 
less time to affect the aircraft position. The reduction in position uncer- 
tainty is also influenced by the higher angular rate of change of the s^ant 
range vector from the VORTAC/TACAN stations to the aircraft. When the 
curves in figure 9 are compared to those in figure 6, it is seen that there is 
an upward shift in figure 9 by a factor of about 2.25. One might expect this 
factor to be 3.0, since the aircraft reference velocity was increased by a 
factor of 3.0 and the Doppler radar velocity measurement error model is a 
linear function of the aircraft reference velocity. This improved performance 
of the velocity estimation is directly related to the improved position esti- 
mation which was explained above. 

Effect of course corrections on the rms velocity estimation error— The 
velocity estimation error for an attitude uncertainty of 0.50“ and 0.25“ is 

shown in figure 9. These two curves each show a peak at the times of both the 

first and second maneuvers of each course correction. In general, the first 

peak is more pronounced than in figure 6, even though a Doppler radar data-set 

sample is processed by the Kalman filter immediately following the first ma- 
neuver. The smaller reduction is a considerably different result than that 
obtained in the subsonic case because the velocity uncertainty in the Doppler 
radar along-track component was three times larger than for the subsonic air- 
craft and will, therefore, have less weight in forming the new state estimate 
in the Kalman filter computations. 

Individual case estimation performance— The total estimation error, r, 
for an individual case is shown for the two attitude uncertainties. As ex- 
pected, most of the minimum points of r occur at VORTAC/TACAN data-set samples. 
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Control of ajtoaZ slate- The total -ms position estimation error 
. IS compared in figure 10 to the total actual rms position 0001^0^’ 

O Total position error (rrms) from the great circle reference course 
standard deviation (la) 

P Along track component of f^ms 

A Cross-track component of frms 

O Total error in the estimation of position ( ) with respect 

to the actual position 

Data sampled: Vorfac/Tacon,doppler rador, barometric altitude 
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also compares ^ 0.50“ and 0.2S“. This figure allows a compar- 

ison of the results for the two attitude uncertainties and a comparison of r 
and for each attitude uncertainty. These results show that the rms de- 

viation of the actual path from the reference path varied from about 3.1 km 
to about 12.1 km for the 0.50“ attitude uncertainty, and from about 2.1 km to 
about 7.4 km for the 0.25“ attitude uncertainty. 


En route timtng errors— The general en route timing performance can be 
seen from figures 9, 10, and 11. Using the en route timing scale, one can 
find the total rms time error, the along-track and cross-track time error, 
and the total uncertainty in the time error at any time along the reference 
path. In general, there would be no difficulty in adhering to fixed time 
schedules along the entire flight path to within a few seconds. On this 
supersonic San Francisco to New York flight, the approximate range of the rms 
en route timing error, , and the rme uncertainty in the knowledge of this 

error, are given b^ow. 


Attitude uncertainty. 

sec 

A./ 

fit , sec 

deg 

rms 

rms* 

0.50 

3.7 - 14.5 

0.6 - 1.6 

0.25 

3.7 - 8.9 

0.4 - 1.1 


This table shows, as in the subsonic case, that the major effect of increas- 
ing the attitude uncertainty is to increase the range of fit 

rms" 

If it is desired to reduce the en route timing error to a minimum at 
some point along the reference path, a course correction can be added at the 
proper time in advance of this point so that r^^g will be a minimum at the 
desired time. Figure 11 shows a flight where a course correction was added 
at 74 min to ensure that the en route timing error at 80 min would be small. 
In this figure, r, the individual case total position error, and ^^ms 
plotted showing the 0.50“ attitude uncertainty situation; curves are plotted 
both with (circled data points) and without the added correction, so that a 
comparison can be made between the results obtained. This comparison illus- 
trates some of the versatility of the system. Not only have the en route 
timing errors at any point along the reference path been maintained at a 
reasonable level, but these errors can be reduced to a minimum at any desired 
time y adding a new correction or by adjusting the time of execution of an 
existing course correction. 


New York City to London - Subsonic Cruise 

In this case, the aircraft is cruising at a reference velocity of 
1,000 km/hr at an altitude of 10.668 km (35,000 ft). The data-sets processed 
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by the Kalman filter were from Loran-C, Doppler radar, and a radio altimeter. 



Figure 12.— Comparison of r , p and for A.U. = 0.25 and 0.50®. 

® ^ rms vms 
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initial Kalman filter settling period, the approximate range of the rms posi- 
tion estimation error, r , and the rms velocity estimation error, v is 
shown in the table belowV 


Attitude uncertainty 

^ 1 
r , km 

V , km/sec 

deg 

rms 

rms 

0.50 

0.30 - 2.2 

8.5 - 13.0 

0.25 

0.15 - 1.6 

5.0 - 9.0 


The table shows that uncertainty in the attitude reference has a more pro- 
nounced effect on velocity uncertainty than on position uncertainty. 

Effect of Loran-C data-set samples— The accuracy of the Loran-C position 
fixes is dependent upon the relative locations of the aircraft and the three- 
station transmitting triad. When the accuracy of the Loran-C position fixes 
is ^oor, the aircraft is said to be in an area of poor coverage. In the plot 
^rms figure 12, there is a large hump during the period from about 50 
min to about 110 min time en route. This hump occurs because the aircraft is 
in an area of poor coverage during this period of time. The poor coverage 
does not have any noticeable effect on v^g because velocity estimation 
depends mainly on the Doppler radar. An on-board system could make some 
improvement in estimation performance in the poor coverage area by modifying 

the sampling schedule to increase the number of Loran-C data-sets in this 
region. ^ 

oouvse Qovreations on the rms veloeity estimation error- The 
effect of course correction execution errors are seen at both the first and 
^econd maneuvers of each course correction for the 0.50® case. The peak in 
second maneuver is generally reduced from that of the first 

the vertical velocity uncertainty component is reduced by 
ter of altitude and altitude rate information by the Kalman fil- 

ter immediately following the second maneuver. For the 0.25® case the 

opposite effect is seen in the third and 

the apparently caused by the fact that 

the vertical velocity uncertainty is always less for this case and the effect 

in ?Sroverafl'ie}n"> " relatively small impro^emLJ 

velocity uncertainty. For the third and seventh course cor- 

ections, v^g was dominated at the first maneuver by good Loran-C fixes. 

error’^'^^''if!^?L^rr performance- The individual case estimation 

error’ of’ abouif 2 ^ Loran-C coverage shows a peak estimation 

^ unusually large throughout the time period 

iHhfarL o/nir ''' performance in position es?imatL 

PYnp^' 1 -dri coverage is consistent with the error which would be 

the timl uncertainty in the position fixes which existed at 
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Control of the aatual states- The rms position estimation for an attitude 
uncertainty of 0.50*’ is compared in figure 13 to the total actual position 
deviation from the reference position, along- track and cross- 

track components of rxms- This figure, as expected, shows that estimation is 
better than control . 

O Total position error (frms) from the great circle reference course, 
standard deviation (lo*) 

□ Along track component of rrms 
A Cross-track component of Trms 

O Total error in the eslimation of position ( Tfirs) with respect 
to the actuol position 

Data sampled: Loran C, doppler radar, radio ottitude 



Figure 13.— rms errors in position and position estimation, subsonic 
case (1000 km/hr). New York City to London. A.U. = 0.5°. 

Present standards for the North Atlantic corridor for lateral deviation 
from the route centerline is 111.1 km (60 n. mi.), except for reductions to 
18.5 km (10 n. mi.) when within 211.1 km (114 n. mi.) of a navigation aid. 

En route timing (along-track) separation standards are apparently ill-defined, 
but are on the order of 20 rain between aircraft (ref. 9) . 

Figure 13 shows that the seven course correction schedule keeps the act- 
ual path (ry^jjg) within 12 km, which is well within the most stringent distance 
standard of 18.5 km (10 n. mi.). 

The flight of the aircraft through an area of poor Loran-C coverage was 
known to occur and previous results had shown that with evenly spaced course 
corrections in this area, the course control could be expected to be poorer 
than when made outside this area. In fact, the poorer course correction per- 
fonricince was observed to result in occasional large actual deviations from 
the reference path in excess of 27 km. This is not an excessive deviation. 
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Nevertheless, it was decided to investigate the effect of delaying the second 
course correction on the overall performance. Figure 14, which compares v 

^ vms 


Subsonic aircraft 
New York City to London 

0.50 deg 

0.25 deg 




Figure 14.- Comparison of and r for A.U. = 0.25 and 0.50". 
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and r for the two attitude uncertainties, shows the results of delaying the 
second course correction to avoid a correction in the area where position 
uncertainty was about three times as large as in adjoining areas. This 
ure shows that although avoiding corrections during a period of larger position 
uncertainty did keep the system from making poor corrections, there was a pen- 
alty due to the relatively long spacing between the corrections in the poor 
coverage area as compared to the spacing in the other areas. This penalty, 
as shown by both and r, is a relatively large deviation of the actual 

position from the reference position, but still results in better performance 
than that observed when corrections were made with a relatively high position 
uncertainty. 

At 246 min, the first maneuver of the sixth course correction is initi- 
ated and illustrates what can happen to r Cfoi" the 0,5® runs attitude error 
individual case) when there is an unfortunate combination of moderate estima- 
tion error and relatively large course correction execution errors. The 
result in this case was that the first maneuver implemented a course correc- 
tion which caused the actual path to cross the reference path, resulting in 
an increase in the actual position error from 4 km to 9.6 km before the 
second maneuver began to return the actual path back to the reference path. 

The rms deviation of the actual path from the reference path varied from 
3.4 km to 11.6 km for an attitude uncertainty of 0.50® and from 2,4 km to 
9.6 km for the 0.25° case. 

En route timing errors— The en route timing errors can be read from the 
time scales on figures 12, 13, and 14. The individual case in figure 14 
shows that due to the poor Loran-C coverage, there was a maximum en route 
timing error, of about 54 sec. The approximate range of the rms en route 
timing error, uncertainty in the timing error, for 

the two attitude uncertainties are given in the table below. 


Attitude uncertainty. 

, sec 

, sec 

deg 

rms 

rms 

0.50 

12.2 - 41.8 

1.3 - 8.0 

0.25 

8,6 - 34.6 

0.5 - 5.8 


As in previous tables, it is clear that en route timing errors are better 
estimated than controlled. Also, while there is a general reduction of 
and as a result of reducing the attitude uncertainty, the magnitude of 

the reductions in comparison to the range indicates the effect is relatively 
minor. The table shows very good en route timing performance when compared 
to the en route timing separation standards, which are on the order of 20 min* 
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The scheduling of a course correction to minimize the error at the ref- 
erence destination time of 300 min is shown in figure 14. This process 
results in an error of about 8 sec for the attitude uncertainty of 0.50* in 
the individual case. 


New York City to London - Supersonic Cruise 

For this case, the aircraft is cruising at a reference velocity of 3000 
km/hr at an altitude of 21.336 km (70,000 ft). The data-set samples processed 
by the Kalman filter were from Loran-C, Doppler radar, and a radio altimeter. 

State estimation perfomanoe— The effect of the two attitude uncertain- 
ties on the rms position, and velocity, estimation error and on 

the position estimation error for an individual case is shown in figure 15. 

The approximate range of these variables for the two attitude uncertain- 
ties are given in the table below. 


Attitude uncertainty, 
deg 

1 

r , km 

ifTiiS 

'O .., . » km 
rms ' 

0.50 

0.2 - 1,5 

18.0 - 29.0 

0.25 

0.13 - 1.2 

6.3 - 19.0 


This table shows that the effect of attitude uncertainty is quite pro- 
nounced on but that there is little effect on In addition, the 

range of Uj^g is much greater than in the subsonic case because of the greater 
speed of the supersonic aircraft. 

Effect of Lovan-C data-set samptes— In figure 15, the hump due to poor 
Loran-C coverage, which was so evident in the subsonic case, is somewhat 
reduced for the supersonic aircraft. This behavior results from the time 
periods between data-set samples being 1 min instead of 3 min in the subsonic 
case, so that errors have less time to accumulate. 

Effect of couvse co2*nections on the rms veXocity estimation otvot — The 
effect of course correction execution errors is seen at both the first and 
second maneuvers of each course correction. The peaks in at each maneu- 

ver are much more evident (particularly when the attitude uncertainty is 0.25*) 
than for the subsonic aircraft shown in figure 12. This effect is due to the 
sensitivity of the course correction execution errors to the aircraft velocity. 

Indiv%dual case estimation performance- The individual case position 
estimation error, v, during the period of poor Loran-C coverage shows a sit- 
uation where the 0.25 attitude uncertainty produced estimates which were 
worse in some instances than those for the 0.50* case (fig. 15). In individual 
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cases, these events are to be expected and, in general, can be attributed to a 
series of unfavorable random numbers being generated by the simulation program 
and to the Kalman filter weighting matrix giving more weight to these errors 
in the 0.25“ attitude uncertainty case than in the 0.50“ attitude uncertainty 
case, as would be expected. 

Contvol of the actual state— The ims position estimation error for an 
attitude uncertainty of 0.50“ is compared in figure 16 with the total position 

O Total position error (frms) from the great circle reference course, 

Standard deviation (lo-) 

□ Along track component of Trms 

^ Cross -track component of 

O Total error In the estimation of position (7rms) with respect 
to the actual position 

Data sampled: Loran C, doppler radar^ radio altitude 



Time en route, min 

Figure 16.— rms errors in position and position estimation, supersonic 
case C3000 kra/hr) . New York City to London A.U. = .5“. 


deviation from the reference, and the along-track and cross-track com- 
ponents or This figure shows that the execution of only three course 

corrections is sufficient to keep the aircraft position errors within 18.5 km 
(.10 n. mi.) (1 a). In the along-track direction, the rms position error is 

expected to be less than about 2.9 km and in the cross-track direction, less 
than 12,5 km. 


Comparison of figure 16 with figure 13 shows that for the supersonic 
case, the along-track position control is much better than the cross-track 
control but that, in the subsonic case, the control performance was 
about the s^e for both. The reasoning which explains the observed perform- 
ance IS as follows: the buildup in position error between course corrections 

is due principally to velocity errors following the previous course correc- 
tion and, to the first order. 
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= (Ar, 7 ,) >< rms speed control error 
'■ AT'’ 2ms ^ kjk-l’ ^ 

and 

^^^CT^rms " k-l^ heading control error 

where and AJ^^y are the 2 ms along-track and cross-track errors, V is the 

aircraft speed, and LT]^ j^_i is the time interval between the second maneuver 
of correction of /c -1 anS the first maneuver of correction k> 

For both subsonic and supersonic cases, the 2 ms speed control error 
(5 km/hr) and heading control error (attitude uncertainty) are assumed to be 
the same. Thus the buildup in along-track error is due solely to the course 
correction scheduling while the cross-track error buildup is due to 

the contribution of groundspeed and the correction schedule, Vi^T-j^ or 

equivalently, the distance between corrections. 

For the supersonic case, Ar^ was generally smaller than for the sub- 
sonic case; for this reason, better along-track control performance is obtained 
for this case. However, the distance between corrections is slightly longer 
for the supersonic case, with the result that the cross-track control perfor- 
mance is slightly poorer. This effect can be seen by comparing figures 13 and 
16. 


In figure 17, and r are compared to show the effects of the two atti- 
tude uncertainties. This figure also allows the comparison of with r for 

each of the attitude uncertainties considered. The upper curve shows r for an 
individual case; the second maneuver of the second course correction (0.25° 
attitude uncertainty) had such large execution errors that the actual devia- 
tion from the reference reached 26.4 km before the first maneuver of the third 
course correction was executed to return the actual path to the reference path. 
This correction was necessarily large and implied another large correction at 
the time of the second maneuver. If a larger number of course corrections had 
been scheduled, the chance of such large deviations from the reference would 
have been reduced. In general, the rms deviation of the actual path from the 
reference path varied from about 2.9 km to about 12.7 km for the 0.50° atti- 
tude uncertainty and from about 1.9 km to about 7.6 km for the 0.25° attitude 
uncertainty. 

En route timing errors— The en route timing errors can be read from the 
time scales on figures 15, 16, and 17. Because of the large course correction 
execution errors at the second maneuver of the second correction, the indi- 
vidual case en route time error in figure 17 reached a maximum of about 32 sec. 

The approximate range of the 2 ms en route timing error, and the 

rms uncertainty in the timing error, 6 * 2 Ws» attitude uncertainties are 

given in the table below 


Attitude uncertainty. 

At , sec 

rms 

6 t , sec 
rms 

deg 

0.50 

3.6 - 15.1 

0.20 - 1.7 

0.25 

2.3 - 9.1 

0.15 - 1,4 
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As in the subsonic case, the table shows the en route timing errors are 
better estimated than controlled. In addition, the table shows a substantial 
improvement in the time control, Strma* when the attitude uncertainty is 
reduced. Comparison of with the subsonic case shows that improved con- 

trol is achieved due to the better along-track position control in the super- 
sonic case. 


Replacing VORTAC/TACAN With Loran-C 

Loran-C was available during about the last 25 percent of the San Fran- 
cisco to New York City flights. In figure 18, the effect of replacing VORTAC/ 
TACAN data-set samples starting at 64 min is shown. The quantity is 



Figure 18.— Comparison of rma estimation error using Loran-C 
and VORTAC/TACAN after 64 min. 


compared both before and after the change to Loran-C. This comparison shows 
a definite reduction in occurs after 64 min time en route, thus demon- 

strating the superiority of the Loran-C position fixes. 
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Effect of Increased Data-Set Sample Rate Before 
Course Corrections 


A short investigation was made to determine what improvement, if any, to 
system performance could be achieved by reducing the position estimation error 
just prior to the execution of a course correction. Two methods were used on 
the supersonic San Francisco to New York City flight. The two methods were: 

(1) retaining the normal sequence of data-set types but doubling the sample 
rate just prior to the first maneuver and retaining this rate until after the 
second maneuver for each course correction; and (2) using the same strategy as 
in (1), but making every other data-set a VORTAC/TACAN data-set and alternat- 
ing the remaining data-sets between Doppler radar and the radio altimeter. 

The results obtained showed that neither of these two methods had more than a 
negligible effect on overall system performance {r^g) , even though the posi- 
tion estimation error was slightly reduced. This effect was observed 

because the limiting factor in system performance is the course control achiev- 
able by the fixed-course-correction-schedule strategy. The dominance of con- 
trol errors over the effect of estimation error on overall system performance 
means an improvement in position estimation will have no more than a negligi- 
ble effect on overall statistical performance. Some individual cases could be 
expected to show improvement, as was seen in the case where a course correc- 


tion was made when the aircraft was in an area of poor Loran-C coverage. The 
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auvctiiLiigfcj ur Dex:i:er szaze esrimation to improve 
system performance would be to improve the course control. One method, 
described earlier, for improving performance would be an adaptive system which 
makes corrections at every opportunity, providing the course correction thres- 
hold is satisfied. 


The effect of reduced velocity estimation error was not investigated 
because early investigations into data-set sample mixing had shown that it 
could not be improved more than a negligible amount by: (1) increasing the 

Doppler radar measurement frequency; or (2) keeping the number of data-sets 
the same and increasing the ratio of Doppler radar measurements to the other 
measurements. 


Effect of the Correlation Decay Coefficient, a 

The correlation decay coefficient, a, in the system dynamics model (appen- 
dix C) was chosen to be 60/hr in all the data presented so far. To demonstrate 
that this choice was not critical, the New York City to London supersonic 
flight was repeated with a = 30 and 120/hr. Since it was assumed that 
would show the most change resulting from different choices of a, this quan- 
tity was compared for a - 30, 60, and 120 in figure 19, The figure shows that 
the results are not sensitive to the choice of a. 


Summary of Course Correction Changes 

The individual case commanded changes in both speed and direction; these 
are given in table 5 for each maneuver of each course correction. The 
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0 10 20 30 40 50 60 70 80 

Time, min 


Figure 19.— Comparison of rms position estimation error for a = 30, 60, 120. 

corresponding statistical changes are also presented for comparison. The sum- 
mation of magnitudes and the algebraic summation of changes in velocity mag- 
nitude and direction are also given for each attitude uncertainty. 

When the statistical entries in table 5 are compared for attitude uncer- 
tainties of 0.50“ and 0.25“ for a particular maneuver, it is seen that for 
0.25“: (1) the changes in magnitude are only slightly less; and (2) the 

changes in direction are of the order of 0.6 times that obtained for the 0.50“ 
case. 

The speed corrections are nearly the same at each maneuver for the two 
attitude uncertainties because they depend, to a first order approximation, 
only on along-track position errors and are independent of estimation and atti- 
tude (heading control} uncertainties, as was shown in the "Course Corrections" 
section of the text. It was also shown in the same section that heading cor- 
rections are, to a first order approximation, only a function of the cross- 
track error. Thus, it is seen that the major influence of heading (attitude) 
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TABLE 5.- COURSE CORRECTION CHANGES 


(a) New York City To London 


o 

■N 

4-J 

O 

o 

4? 

o 

u 

Q) 

12 

3 

(3 

Ar 

0> 

3 

<D 

# 

Supersoni< 

: aircraft 

Subsonic aircraft 

Command ec 
corrc 
(individL 

velocity 
ictions 
lal case) 

ms commanded velocity 
corrections 

Commanded velocity 
corrections 
(individual case) 

ms commanded velocity 
corrections 

A.U.^ 

= 0.5 

A.U. = 

= 0,25 

A.U. = 0.5 

A.U. = 0.25 

A.U. = 0.5 


A.U. - 0.5 

A.U. - 0.25 



As 

Aij; 

As 


As 

Aij^ 

As 


As 

mm 

As 

Aijj 

As 

Ai); 

1 

1st 

2nd 

-12.50 

11.18 

1.49 

-1.00 

-12.08 

12.59 

1.09 

.00 

39.92 

32.06 



1.21 

1.00 

-11.24 

16.82 

0.43 

0.53 

-11.38 

20,12 


27.44 

18.51 

0.62 

.42 

27.38 

18.44 

0.51 

2 

1st 

2nd 

-7.92 

12.24 

1.74 

-2.22 

-19.39 

27.07 

-2.21 

3.23 


IQIQ 

■AM 

■allEl 

mBm 

.00 

.00 
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jjHQQn 

.16 

.18 

3 

1st 

2nd 

38.95 

-23.82 

-2.64 

3.73 

41.37 

-42.11 

-7.17 

6.29 

41.12 

34,95 
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2.35 
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.66 

mKjkKkm 

mSWSM 
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.44 
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1st 
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■Mg 
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.46 



5 
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.00 

12.50 



m 

■ullSI 

.48 

.43 


mm 

6 

1st 

2nd 









.00 

-5.81 



BJUJ 

Kil 

.50 

.46 



7 

1st 

2nd 










■IRl 


.00 

-.76 




IH 

summation ox 
magnitudes 
a. 

106.41 

12.82 

154.61 

19.99 

220.72 

15.41 

207.42 

8.42 

97.68 

25.33 

277.69 

— 

17.77 

335.90 

6.87 

355.62 

4.63 


= me Attitude (heading control) uncertainty, deg. 
= Change in ground speed, km/ hr. 

= Change in heading, deg. 

















TABLE 5.- COURSE CORRECTION CHANGES - Concluded 


Cb) San Francisco To New York City 



■ 



Supersonic aircraft 





Subsonic 

aircraft 




■M 

o 

■ 

Commanded velocity 

rms 

commanded velocity 

Commanded velocity 

rms commanded velocity 




corrections 



corrections 



corrections 



corrections 


o 

o 

H| 

(individual case) 






(individual case) 






Qi 

Hi 


IB 





A.U. = 

0.25 

A.U. 

= 0.5 

A.U. = 

0.25 

A.U. = 

0.5 

A.U. = 

0.25 

cl 

H 



hs 



A4) 

As 


As 

AV; 

As 

A\[j 

As 

A()j 

As 



la 

HHH 

■IM 





34.57 

1.13 

-10.63 

3.08 

-11.52 

2.44 

21.89 

0.48 

21,87 

0.07 

1 

HBH 

iHRnl 

kH 


■in 



26.11 

.96 

6.29 

.00 

5.27 

.50 

14.88 

.31 

14.86 

.25 


1st 

14.41 

.00 

7.60 

-4.46 

43.26 

2.63 

40.04 

1.60 

36.77 

-4.72 

27.61 

-1.55 

28.91 

.49 

28,67 

.40 

2 

2nd 

.00 

.47 

-13.52 

2.86 

37.07 

2.25 

33.49 

1.40 

-36.64 

2.22 

-24.38 

-1.44 

21.62 

.43 

21.40 

.34 


1st 

.00 

-.85 

83.17 

2.54 

36.29 

2.34 

33.86 

1.42 

12.54 

4.72 

-17.09 

4.12 

29.14 

.48 

28.98 

.34 

3 

2nd 

-8,08 

.00 

-86.18 

-2.05 

30.75 

2.01 

28.06 

1.25 

-7.97 

-1.22 

18.70 

-3.23 

22.01 

.40 

21.83 

.29 


1st 









12.39 

-1.09 

13.89 

0.25 

28.78 

.46 

28.72 

.31 

4 
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-5.98 

.00 

.00 

-.82 
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.40 
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.00 

2.30 

14.61 

.00 

28.82 

.31 

28.81 

.28 

5 
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-11.83 

-2.18 

-19.91 

.25 

22.75 

.38 

22.66 

.25 

Summat 

ion of 

28,32 

3.90 

190.47 

13.08 

212.04 

12.90 

196.13 

7.76 

141.04 

21.58 

152.98 

14.60 

240.77 

4.14 

239.67 

3.10 
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^A.U. = rme Attitude (heading control) uncertainty, deg. 
= Change in groundspeed, km/hr. 

= Change in direction (heading), deg. 







uncertainty is on the cross- track position control; changes in heading uncer- 
tainty result in nearly linear corresponding changes in the cross -track com- 
ponents of and u 

^ rms rms 

The commanded changes (individual case) are signed values to indicate the 
direction of the change. These commanded changes exhibit some randomness in 
magnitude when the first and second maneuvers of each course correction are 
compared. The signs of the heading and speed commands for the second maneuver 
will almost always be opposite to those of the first maneuver. This effect is 
observed because, ideally, the first maneuver of the course correction causes 
the aircraft to intercept a nearby point on the reference trajectory, matching 
both position and time at that point. The second maneuver then matches the 
reference velocity vector at the intercept point. This strategy requires rel- 
atively large changes at the first maneuver and makes a second maneuver manda- 
tory in order to remove most of the speed and heading changes applied at the 
first maneuver. The advantage of this strategy is that timing relative to a 
reference path is controlled at any time along the reference path. This en 
route time control advantage is obtained by using speed (thrust) changes in 
the correction maneuvers in addition to the heading changes necessary to con- 
trol the cross-track position errors. 


■+■ 1-1 v-v -m 1 mV. n-i/x ^ -C ..wvw. „ l_ - 


tiler vy a. nxea-scneauie 
or by use of an adaptive control system, will improve the en route time and 
position performance of the system. When this process is accomplished, it can 
be expected that there will be an increase in the number of individual (speed 
or heading) commanded changes which will be so small that they will be below 
the threshold of the course correction system and will be skipped. In effect, 
it is these speed and heading thresholds in the automatic control system which 
establish the lower limit of position and en route time control which can be 
achieved. In this report, individual commanded changes which were below their 
respective thresholds and were skipped by the automatic control system are 
shown as zeros in table 5. 


Each column in table 5 has a summation of the magnitudes of the values 
appearing in that column. This summation is an indication of the total throt- 
tle activity required for time control. The required thrust changes for each 
maneuver are proportional to As. It can be seen in table 5 that, for each 
case, the rms speed changes are approximately the same for every speed correc- 
tion; hence, rms thmist changes are about the same each time the speed is 
corrected. 


Comparison of the results in table 5 for the two attitude uncertainties 
indicates that the total rms As requirements are independent of attitude uncer- 
tainty. It can also be seen from the individual cases that the speed change 
in the second maneuver approximately cancels the speed change of the first 
maneuver at each course correction. The cancelled portion is due to the time 
correction and the residual speed change is due to correction of the initial 
error in speed from the reference speed. 

The summation of the magnitudes of the heading changes indicates the 
total heading control activity used to control the cross-track position. The 
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attitude uncertainty has a pronounced effect on the heading control activity, 
as expected. 


RESULTS AND CONCLUSIONS 


The results of this study show that a navigation system such as the one 
described in this report has excellent en route time control capability 
throughout the entire flight of subsonic or supersonic aircraft traveling 
along either of the great-circle routes considered. It is shown that good 
system performance throughout the entire flight requires both good state esti 
mation and good control over the actual state. 

The results summarized below are taken from the digital simulation after 
sufficient time has been allowed for the system to "settle out." All summary 
data presented in this section will be for an attitude uncertainty of 0.50° 


Summary of Results 

The basic assumptions which were made for the transcontinental flight are 
summarized in the table below 



Supersonic 

Subsonic 

Altitude 

21.336 km (70,000 ft) 

10.668 km (3.5,000 ft) 

Speed 

3,000 kra/hr 

1,000 km/hr 

Course correction 

Every 30 min (3 course 

Every 45 min (5 course 

schedule 
(see table 1) 

corrections) 

corrections) 

Data types 

TACAN/VORTAC, Doppler, radio, 
barometric altimeter 


Data sampling 

Sequential, once every min 

Sequential, once every 
3 min 


Similarly, the assumptions vv^hich were made for the trans-Atlantic flight are 
summarized in the following table. 



Supersonic 

Subsonic 

Altitude 

21.336 km (70,000 ft) 

10.668 km (35,000 ft) 

Speed 

3,000 km/hr 

1,000 km/hr 

Course correction 

Variable — approximately 

Variable — approxi- 

schedule 

every 25 min (3 correc- 

mately every 45 min 

(see table 1) 
Data types 

tions) 

Loran-C, doppler radar, 
radio altimeter 

(7 corrections) 

Data sampling 

Sequential, once every min 

Sequential, once every 
3 min 
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The results obtained for the transcontinental and the trans-Atlantic flights 
are presented in terras of an approximate range for the following variables: 

(a) rms en route timing error, uncertainty in the en route 

timing error, 6ii;>ms> deviation from the reference path, (d) rms 

position estimation error, (e) rms commanded speed correction. These 

variables are tabulated below for an attitude uncertainty of 0,50° 



St , sec 
rms 

St , sec 

rms 

r , km 

rms 

r , km 

ms 

As 

ms 

Transcontinental flight 

Supersonic 

3.7 - 14.5 

0.6 - 1,6 

3.1 - 12.1 

0.35 - 1.25 

28,4 - 43.3 

Subsonic 

13.5 - 30.0 

1.8 - 6.0 

3.8 - 8.3 

0.50 - 1.65 

14.9 - 29’. 1 

Trans -Atlantic flight 

Supersonic 

3.6 - 16.4 

0.20- 1.7 

2.9 - 12.7 

0.20 - 1.50 

32.1 - 41.1 

Subsonic 

12.2 - 41.8 

1.3 - 8.0 

3.4 - 11.6 

0.15 - 2.30 

13.8 - 42.5 


Conclusions 

1 . The maximum rms actual position deviation from the reference position 
is dependent upon the number and spacing of the course corrections. 

2. The system performance is limited by the execution errors in the con- 
trol system and external forces in the physical environment. Thus, the real- 
izable minimum values of r^^Q at the intercept point on the reference path is 
significantly larger than r^^s the aircraft trajectory is better estimated 
than it can be controlled. 


3. The minimum ms value of can be reduced slightly by increasing 

the data-set sample rate just prior to making the course corrections. 

4. A minimum rms actual position deviation from the reference (and rms 
actual en route timing error) may be moved to any desired point along the ref- 
erence flight path by changing the course correction schedule so that a mini- 
mum rms deviation occurs at the desired time. 

5. Loran-C position fixes are more accurate than those obtained from 
VORTAC/TACAN position fixes. 

6. The system performance is relatively insensitive to the choice of 
value used for the correlation decay coefficient, a, in the dynamic system 
model within the range from 30 to 120/hr. 

7. The supersonic aircraft guidance and navigation system in this report 
could meet the lateral separation standards presently specified by the FAA 
over the continental United States by an increase in the number of course cor- 
rections from three to four or five. The en route timing error standard was 

easily met. Thirty sec or less was obtained with the schedule used in this 
study. 

8. Along a great-circle route in the North Atlantic, the performance of 
the guidance and navigation system in both the lateral and longitudinal 
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directions was such that the lateral position error was less than 13 km and 
the en route timing error was less than 42 sec on all flights. 

9. It is postulated that a modification to the course control strategy 
would provide an adaptive control system capable of substantially improved 
performance. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. , 94035, May 10, 1973 
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APPENDIX A 


EFFECT OF GEOMETRY ON LORAN-C POSITION UNCERTAINTY 


To obtain a position fix, a Loran-C receiver is used to obtain two simul- 
taneous time difference readings from a master and two slave stations in a 
particular group or chain of stations identified by a basic and a specific 
pulse repetition rate. The time difference readings define two hyperbolic 
lines- of-position (LOP's) whose point of intersection is the desired position 
fix. The accuracy of this fix is dependent upon the geometric location of 
the receiver with respect to the master and two slave stations and on the 
receiver itself. The geometric location of the stations used in this report 
is that of January 17, 1969 (reference 5), and is shown in figure 2 in the 
text. When a particular station group is chosen, the effect on position 
accuracy due to that particular geometry with respect to the receiver is 
fixed and arises because of the hyperbolic nature of the LOP's and the fact 
that the two LOP's do not, in general, cross at right angles. 


The errors in measuring the two time differences in the time of arrival 
of the pulses from the selected master and two slaves, say x and y, are 

assumed to have Gaussian distributions with variances o ^ and c 

X y 

The position fix is desired in terms of latitude and longitude and it is 
therefore necessary to find the components of cr^,^ and in these spherical 
coordinates. The method for finding these variances is" given in reference 10 
but was modified somewhat to obtain the results given in reference 3 for an 
intermediate orthogonal coordinate system (V,Z). The results in terms of 
equivalent distance units are: 


2 _ 


C-2 


w 


4 sin <5 
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, O 

sin^ n 
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X 

4 sin^ n 
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2p a a cos 6 

cos2 5 + fy ^ y _ + 
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o cos 6 P o 
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sin^ n 


cos^ 6 + 


X 


sin^ n 


y 


sin sin 


1/2 


(A3) 


where C is the velocity of light in kilometers/microsecond and and are 
one-half of the angles between a line drawn from the point of LOP intersection 
to the master station and a line from the same point to slave station x and 
to slave station y, respectively. (See sketch (Aa)0 The LOP crossing angle, 
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<5, is given by 



Sketch (Aa) 


The correlation coefficient, is due to both time differences being mea- 

sured with respect to the same master station pulse. Reference 10 gives 0.3 
as the usual value chosen for p^. 

The results in equations (Al), (A2), and (A3) can be expressed as a co- 
variance matrix, Q, which represents an error ellipsoid in the 2 coordinate 
system. Thus 



(AS) 


Essentially following references 3 and 10, the major axis of the error ellip- 
soid lies along a line from the position fix (LOP crossing point) to the 
master station, with its center located at the position fix. Let two unit 
vectors, £/f/ and i/g, lie along the major axis and the minor axis, respectively, 
of the error ellipsoid as shown in sketch Aa. This sketch shows the relation- 
ship of the error ellipsoid to the station geometry and includes: 
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(1) two unit vectors, % and U^, which are directed to the north and east, 
respectively, from the position fix; (2) typical directions to the master and 
two slave stations; and (3) th^ angles na; and wy, which in equation (A4) 
define the LOP crossing angle, 6. From this sketch, the angle, c, between 
the W, Z coordinate system and latitude, longitude coordinate system ri/p, Un) 
is found from 


cos c = = V^-V^ 


(A6) 


This angle allows transformation of the error ellipsoid into the latitude-, 
longitude coordinate system using a transformation matrix of the form 


h = 


cos t, -sin C 


[sin K, cos 


(A7) 


The covariance matrix of the position fix uncertainty in terms of 
distance squared in the latitude, longitude coordinate system XJ is given 
by: 


Given the magnitude of the radial distance from the center of the earth 
to the aircraft, |i?| , and an estimate of the latitude, ^ transformation 
matrix A is formed which transforms into angular units (radians! squared. 
Thus, 


where 






fsec 4> 


o 


0 

1 


The matrix is used as the measurement uncertainty matrix for Loran-C 
in the Kalman filter and in the generation of the random errors for simulating 
the output of the Loran-C coordinate converter; the output of this converter 
is latitude and longitude, given two time differences from a master and two 
slave stations. 
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APPENDIX B 


TRANSISTION MATRIX FOR THE EQUILIBRIUM FLIGHT PATH 


A general transition matrix will be derived which allows 

the aircraft state at time t = * "z to he found from the aircraft state at 

t - tQ along an equilibrium path. The general form of this transition matrix 
will be such that it will be applicable to all three aircraft states used in 
this report (i.e., the reference, actual, and estimated states), requiring 
only that the transition matrix be evaluated from the appropriate state 
variables. Thus, the computation to find the estimated aircraft state at 
t - tQ ■¥ 1 is simply 

* t; = CBI) 

where 


B(t) 

O 


x(to> = 




In this study, the basic form of equation (Bl) will be used to predict the 
equilibrium path from one data-set sample to the next and for predicting the 
effect of added velocity increments in computing course corrections . The 
prediction is in error due to unknown randomly occurring forces during the 
interval, such as automatic control system errors and wind disturbances. 


The equilibrium kinematics are defined by the assumptions that 

1. The angular momentum of the aircraft with respect to the earth is 
conserved. 

2. Altitude rate is constant over the period t. Equivalently, the air- 
craft is assumed to be flying in a plane (great-circle) with respect to the 
earth, with constant altitude rate, and a groundspeed-altitude relation that 
conserves angular momentum. 


The equi-lihrium pat/z— The motion of the aircraft is described in two co- 
ordinate systems: (1) a right-handed Cartesian system rotating with the 

earth, in which the Z-axis passes through the North Pole and is the axis of 
rotation, the X-axis coincides with the meridian of zero longitude, and the 
Y-axis is at 90* east longitude; and (2) a system consisting of a unit vector, 
which lies along the instantaneous position vector, R(t)\ a unit vector, 
Zj,, which is perpendicular to ij, in the plane of motion and in the direction 

of motion; and a unit vector, i , perpendicular to both i and % , and given 

q ^ p 
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These coordinate systems are shovm in the following sketch (Ba) . 


t>y i 


North Poie 
Z 



The position and velocity of the aircraft may now be defined as 

R ^ A 


V ^ V i * H 
VP 3 " 


(B2) 


Since the altitude rate is assumed to be constant during the period x, then 

v(x) = h(t J = r 
o o 


= r + vji 
o o 


(B3) 


Because the angular momentum 

i?xy = vv I 

p <7 
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is conserved, then 


and 



- v^(yr(y - 


CB4) 


or 


" ^po^o/^^o * 


po 

y\ 4S /V 

The vectors and ipCr) are found as follows: Since = o, the 

aircraft motion remains in the plane defined by Rq and V^, hence iQix) ^ 

and ipo are all in the same plane. Sketch (B-b) shows the re- 
lationship of these two sets of orthogonal unit vectors after rotation in the 
plane of motion through the angle 6. 
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From this sketch, an expression relating the two sets o£ unit vectors may be 
written as follows: 


i^(T) 




cos Q(t) 


sin Q(r) 


-sin Q(t) cos 


ro 


po 


The unknown angle, Q(r), is evaluated from 

e(T) = V^(’i)/r(x) 

which, after integrating, using equations (B3) and (B4) , gives 

V T 

Q(r) = — = riT 

r + r T 
o o 


(BS) 


(B6) 


The state at t = to + t may now be expressed in terms of the initial state 
and the time increment, t. Combining equations (B2) , (B3) , (B5) , and (B6), 
we obtain: 


R(t) 


[- 


w COS nt - — sin t\t 
o r 


1 . i 

Jon 


sin nxt' 


(B7) 


V(x) 


= r pnx cos nx - Ty + n7sin nx 1 fi? + — — sin nx + n cos nx | V 

L J L o o» r Jo 


where 


V • 0 V 

po , '^po 

- — * — , u = — , and n = ^ . — 

o r ’ p + r x 

o po o o o 


In transition matrix form, with Eq and Vo given in terms of the Cartesian 
coordinates shown in the first sketch in this appendix, equation (B7) 
becomes 


r n 

ECxJ 

- ^Ct + X, t ) 

0 o 

ro 

0 

vrxj 
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(B8) 
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where 






COS nt 


I (i sin nT)j 


- — sin riT^J 


|yr|T cos r)T 

[. ('y + n>>sin n-rjl 


1 { i> . 

(0 (j’ 

O \ <3 


(B9) 


Sin ni 


+ n cos ni 


U 


where I is the 3x3 identity matrix. The quantities in i() are functions of 
initial conditions and time and are evaluated for the estimated, actual and 
reference paths using the respective initial conditions for the path in 
question. 
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APPENDIX C 


COMPUTATION OF AND THE MATRICES Ai? AND A? FOR 

WANDER FROM THE PREDICTED PATH 


The Kalman filter is used to process measurement data-sets every At 
hours. Between any two of these measurement data-sets, say at and 

t-j^, the aircraft velocity vector wanders in pitch, heading, and velocity mag- 
nitude from that predicted by the transition matrix for conditions at 
This non-deterministic wander from the predicted path is due to random dis- 
turbances such as changes in the wind, air density, temperature, and so forth. 
The on-board automatic control system is assumed to be capable of limiting the 
excursions due to these disturbances but, because of errors in the measurement 
devices and errors in the automatic control system itself, the aircraft does 
not follow the desired (commanded) flight path. In addition, because of the 
instrument and estimation errors, it was assumed that the on-board system 
would not have a precise measure of deviations from the commanded flight 
path. Thus, only an estimate of the actual flight could be made. Further- 
more, because of the dynamics of the aircraft, its automatic control system, 
and time constants in the measurement devices, the deviations from the 
commanded flight path were modeled as time correlated error sources, as were 
the instrument errors. The statistical model of the deviations of the actual 
state from the reference state will be considered in two parts: (1) non- 

stationary; and (2) stationary. Nonstationary statistics result when a 
course correction is executed. At this time, the Kalman filter processes a 
final data-set and the updated aircraft state is used in the course correction 
computation to find incremental changes in speed, heading, and pitch. The 
automatic control system is directed to change rapidly the aircraft heading 
and pitch by the incremental values and then change the speed by the command- 
ed incremental value. During the short period of time the speed is being 
changed, the uncertainties in heading, pitch, and roll are equal to the short 
period instrument errors of the vertical and heading reference. Since the 
long period errors have a time constant of 60 sec, the error in the control 
system implementation of the incremental changes is assumed negligible. 
Following the completion of the course correction maneuver, the uncertainties 
in the new pitch, roll, heading, and speed increase as the long time constant 
errors come into play. That is, the control system must now maintain total 
pitch, roll, heading, and speed rather than implement an incremental change 
as was done in the course correction. 

The long period errors are assumed to be characterized by the following 
differential equation: 
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CD 






= -a 


+ a 







1 


L _j 


L J 


(Cl) 


where 
in 


re Og, 0 , a , and Og are the long period constants taken from table 2 
the text. A solution to equation (Cl) is 



(1 


-at 


(C2) 


where t is equal to nLt — the number of At intervals occurring since the 
last course correction. When 10 At intervals have passed (n = 10), the 
statistics of the deviations are considered to have exponentially reached 
stationary^ values. Under these conditions, the primed values on the left 
side of equation (C2) are replaced by the unprimed values on the right. 


The error sample generation process used to obtain the actual state 
deviations from the reference state at is essentially four identical 

first order filters driven by individual random noise sources whose variance 
is determined by equation (C2) . The well-known discrete form of the sample 
generation process which simulates the effect of a first order filter is 
given below: 
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;c-l 

s 4 


(C3) 


^This assumption is not strictly true in the case of the heading error 
standard deviation which is actually increasing at 0.1®/hr due to a pseudo- 
gyroscope drift discussed in the text. 
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where the (i - 1, 2, Z, A) are Gaussian, zero mean, unit variance, random 
numbers which are independent of each other and uncorrelated in time. Equa- 
tion (C3) generates error samples which are correlated in time with previous 
samples and is an expression normally used with stationary statistics; that 
is, the variances of the samples generated by this expression are constant 
when the variances given by equation (C2) have reached stationary values. 

In addition, these pitch, heading, and speed samples are measured in an ALT, 
ALONG-TRACK, CROSS-TRACK coordinate system defined by the reference state at 

The roll angle, A(J), has little significance in the simulation results, 
but was carried along for completeness. The limiting action of the automatic 
control system was simulated by exami ning the heading deviation, Aip, given by 
e q^uatio n (C3) and, if larger than ±/ 2,0", was reduced to the nearest 

±/27o^. 

The samples obtained from equation (C3) represent the deviation of the 
actual trajectory from the reference trajectory. The samples represent 
changes in the aircraft trajectory because of variations in such factors as 
wind, air density, control system errors, and so forth. If it were not for 
these deviations, the aircraft would fly along the reference path by virtue 
of the fact that nominal corrections have already been applied to the aircraft 
speed, heading, and pitch attitude by the on-board computer. 


In this study, a vector, AX^, is used to modify the existing actual state 
so that it corresponds to the actual state indicated by the error samples 
obtained from equation (C3) . This procedure is implemented immediately before 
the processing of a data-set by the Kalman filter. 


Derivation of AX^ 

To compute the AX^, used in the text, an expression is required for the 
new actual velocity vector in terms of the error samples which have already 
been found. Thus, the new actual velocity is given by 


where 



(C4) 


(C5) 


is the new groundspeed obtained by adding the speed change to the 

speed predicted by the transition matrix from conditions at 
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Let the predicted actual state at tj^~ be expressed by 


X 


* 


A1 




A2 


(C6) 


*k- 


and let the desired in the same coordinate system be given by 




AX’ 


AX 


Al 


A2 


(C7) 


The required change in velocity can now be written as 




CC8) 


The position deviation from that predicted by the transition matrix from 
conditions at is obtained by assuming that this change is parallel to 

the new velocity vector in the ALT and CROSS-TRACK directions, but equal to 
At times As in the ALONG-TRACK direction. Thus 



AX - U 


As 




(C9) 


All of the elements of equation (C7) have been determined in the reference 
ALT, ALONG-TRACK, CROSS-TRACK coordinate system. To find axa> that re- 
mains is to transform equation (C7) into the X, Y, Z coordinate system. 

Thus the required matrix operation is 




T 0 


0 T 




(CIO) 
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where T is the known 3x3 transformation matrix. 


Derivation of the Matrix Ai? 


Equation (18) in the text shows the addition of a matrix Ai? after the 
transition matrix has been used to update R(ti^.i) to R(t-j^) . The quantity 
Ai? is added to account for the increase in uncertainty in the deviation 
between the actual and reference trajectories following course correction 
maneuvers as given by equation (C2) . Equation (C2) produces nonstationary 
variances following each course correction and reaches a stationary value 
after 10 data-sets have been processed. It is this increase in the variance 
for each At that must be accounted for by Ai?. 


n 


Returning to equation (C2) , the change in the left side as a function of 
is given by 


— 





% 
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(Cll) 


Now, translating into equivalent ALT, ALONG-TRACK, and CROSS-TRACK velocities 
and positions we have 


6o0 I 7^1 At 
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The matrix 


i?y* is given by 




W 


w 


w 


(C13) 


where both and are diagonal matrices and the elements of 

(R^ ) f are the sauare roots of the elements of rJ" . Thus, the diagonal 
elements of (R^ )l/2 are the elements of the 6xi vector 6x*; all other 
elements are zero. 
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To find Ai? it is first necessary to transform into the ref- 

erence ALT, ALONG-TRACK, CROSS-TRACK coordinate system to^obtain R* (t-) . 
Define ^ 


D(t.) and D(t-) 
^ k 


diagonal matrices whose elements are the 
square roots of the diagonal elements of 

R*(t,) and R* (t.) , respectively. 

^ k 


(C14) 


Now, 


hR* = D^(t.) - D^(t~) 
k k 


= D(t.) + R *1/2 _ (CIS) 

k ^ k 

= R* * 2D(t_)R*^/^ 
w ^ w 


The matrix Ai? is found by using the transformation matrix, T, and trans- 
forming AS* as follows: 


T 


AS = 


0 




T 


(C 16 ;) 


Derivation of the Matrix AP 

During the interval ("i^k-l* between measurements, the uncertainty 

in the state estimate obtained at increases because of random changes 

in the wind, air density, and so forth. This increase is assumed to satisfy 
equation (Cl) with the following solution: 




AOg 


»e - 




= 


(l-e ) 

(C17) 

Act 

L 'liJ 


a - \ 




where the a's are the same as those in equation (Bl) and the 
j = 4, 5, and 6, are the square roots of the diagonal elements of the P 
matrix, which has been transformed into the actual ALT, ALONG-TRACK, CROSS- 
TRACK coordinate system. 

In equation (C17) , the column vector on the right is the driving uncer- 
tainty. It is obtained by subtracting the rms uncertainty in the current 
estimate of the aircraft pitch, speed, and heading from the corresponding 
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rms instrument uncertainty. These instrument uncertainties are the values to 
which the pitch, speed, and heading uncertainties would eventually rise if no 
further aircraft state estimates were obtained through the processing of 
navigation information. The matrix Py is the covariance of the increase 
given by equation (C17), Therefore, following the method used to find AP, 
we have 
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from which 


. (1 . 
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Ao Lt 
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0 40^|7j,| 


0 




0 Aa 


0 

0 


0 


Now, using the same approach used to find AP, define 


and DP(t^) = diagonal matrices whose elements are the 
^ square roots of the diagonal elements of 
Pft^j and P* (t-) , respectively. 


(C19) 


(C20) 
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Then, 


AP* = UP'^it.) - 

^ k 

= [DP(t.) * PyV2]2 . VP^(t ) 
K ^ K 

= P + 2DP(t,)P 1/2 
W k ^ 


(C21) 


The required AP is then found by first transforming AP* from the actual 
ALT, ALONG-TRACK, CROSS-TRACK system into the reference ALT, ALONG-TRACK, 
CROSS-TRACK system using the actual angular deviations from the reference 
(found earlier) to compute a transformation matrix IT^. Once this process is 
accomplished, a second transformation is made to the X, Y, Z coordinate sys- 
tem using the transformation matrix T. Thus, 









TT, 
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TTa 
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(C22) 


Derivation of hX„ 

tf 


The vector AJfg, is computed in a manner completely analogous to the 
procedure used to compute AZ^. Therefore, only the differences will be 
emphasized here. 


The angular deviations from the actual path must be found using the 
updated elements of DP(tj^) which are also the square roots of the diagonal 
elements of P* . The deviation angles are then found from 
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where the U^, i = 1, 2, and 3, are Gaussian, zero mean, unit 
random number samples. The estimated quantities are produced 
samples in equation (C2) the deviation obtained from equation 


(C23) 


variance, 
by adding to the 
(C23) . Thus, 
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From this point, the derivation follows that of AX^ exactly, except that 
the estimated quantites found in equation (C24) are used instead of the 
actual quantites found in equation (C2) . When all the elements of (t]^) 

have been found as in equations (C7] through (C9) , then a transformation 
matrix is found from the error angles in equation (C23) . The required 

transformation equivalent to the expression in (CIO) is 


Ax^rt^; = 


TT^ 0 


0 TT^ 


AXg-* (ti^) 


(C25) 
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APPENDIX D 


DERIVATION OF THE MEASUREMENT JACOBIAN MATRICES FOR 
LORAN-C, VORTAC/TACAN, DOPPLER RADAR, AND 
RADIO AND BAROMETRIC ALTIMETERS 


Each time a data-set is processed by the Kalman filter, a measurement 
Jacobian or H matrix is required which relates changes in each measured 
quantity to changes in each of the state variables. These matrices, which 
are evaluated on the reference trajectory, are derived in the following^ 
sections. The basic coordinate system defined by the unit vectors i, j, and 
k is the usual geographic system rotating v/ith the earth, as shown in sketch 
CDa) . 


North Pole 
Z 



LORAN-C 

Loran-C position fixes are obtained by simultaneously measuring two time 
differences between radio signals received from a three-station triad. These 
time differences are then fed to a coordinate converter which, knowing the 
position coordinates of the three-station triad, computes the latitude and 
longitude of the aircraft. The measurement Jacobian matrix, is of the 

form: 
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'^(Lat) d(Lat) 

3 (Long) 
dT~ 


dX 

3 (Long) 

. dx 


3 (Lat) 
dZ 

8 (Long)\ 
3Z 


Each element of is evaluated from: 

Lat = sin“^ ^ 

Y 

Long ~ tan’^ j 
R = + z2;i/2 

Taking partial derivatives, we have: 

-XZd -XZd d} 
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where 


^ R'^d'^ Lj?2j r^x 

d = (X'^ + y2;i/2 
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(D3) 

(D4) 
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VORTAC/TACAN 

For both VORTAC and TACAN, it is assumed that slant range and bearing to 
the ground station from the aircraft are measured simultaneously. 

In sketch (Db), a typical configuration is shown. The slant range 
vector is from the ground station to the aircraft and the bearing angle, 6, 
is measured from true north to the projection of the slant range vector onto 
the tangent plane at the ground station. In practice, bearing is measured 
from magnetic north, but the difference between true north and magnetic north 
is known at each station, and it is therefore assumed that this correction is 
stored in the computer to correct the measured bearing to true north. Also, 
when the on-board equipment measures the magnitude of the slant range vector, 
the on-board computer is assumed to have stored a predetermined bias correc- 
tion which is applied prior to any further computations. 

The following definitions are made for the geometric elements shown in 
this sketch: 

1. The symbol Q is the bearing from the VORTAC/TACAN ground station 
to the airborne receiver, measured clockwise from true north to the vector 

■^1-. * 
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2. The symbol S is a vector from the center of the earth to the 
ground station and is a constant vector for the particular VORTAC/TACAN 
station. 

3. The symbol E is the current vector from the center of the earth to 
the on-board receiver. 

4. The symbol is the slant range vector from the ground station to 

the on-board receiver. It is the vector difference of R and S', that is, 

^ E — S » 

5. The three unit vectors i/„, U^, and Ug constitute a rectangular 
coordinate system. The unit vector Ug, lies along the S vector. The unit 
vector, U^, lies perpendicular to S and in the direction of true north. 

The third unit vector, which is perpendicular to both Ug and U^, com- 
pletes the right-handed orthogonal coordinate system and lies to the east. 
These unit vectors are constant for a particular ground station and are com- 
puted from the following relationships: 
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« = 


s 

!5| 


(D7) 


U = 


s 


(D8) 


U = -; 
n I 


U 

e g 

U >^V I 
e s 


(D9) 


6. The symbol /?□ is the projection of the slant range vector, 5^, 
onto the ground plane defined by the unit vectors K 


^7 = 


1 matrix. 

Hy, is 

of the form: 


38 

38 

38 



dX 

9J 

92 





9|^i?l 


(DIO) 


3Y 

92 _ 




To find the required partial derivatives of the bearing angle, 8, the compo- 
nent of the slant range in the northerly and easterly directions are found 
using the slant range vector and the unit vectors, and U^, as follows: 

N = cos 6 (Dll) 

E = S^-U^ = sin 8 (D12) 

from which 

R 1 . Jj2 + r2 (-[, 13 -j 

n 

The partial derivative of 8 with respect to N is found by taking the 
partial derivative of both sides of equation (Dll) which yields 

If = I (t¥ - i) (‘>1-') 

To evaluate both sides of equation (D13) are differentiated with the 

result 
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CD15) 


9^ i?. 


SO that 


9B 1 / if ^ \ -E 


-sin 6 


(D16) 


Using equations (D12) and (013), the partial derivative o£ g with 
respect to E is found in a manner analogous to that just described. This 
process results in the expression 

96 if cos 6 fn- 


Since 


lies in a plane perpendicular^ to the local vertical at the 


ground stations, it has no component along U . Therefore 

s 


The upper row of is a row matrix^of^the coefficients of the gradient of 

the scalar bearing angle, 6 , in the i, k coordinate system. Since the 

gradient of a scalar is invariant under a unity orthogonal transformation, 

the gradient of 6 in the t, 3 ^ k system may be found^by^transforming the 

gradient of 6 in the U^, system into the i, 3, k system. The 

required transformation matrix is the transpose of a matrix composed of the 

unit column vectors U , U , and U , which have already been found. Thus, 

n e s ^ 


f£6 16 96. 1 ^ ril 96 96 1 r 
[9X 97 9Z J [ 9^7 9£’ aFj I 


U U U 
n e s 


is the required matrix equation for the upper row of Hy. 

To find the elements for the lower row of Hy, the slant range may be 
expressed as 

= |i? - 5| = [(X - Si)^ + (7 - 82)^ + (Z - 53 )^] ([ 

where 


R = Xi yj * Zk 


(D21) 
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and 


S = Sli 

+ S2<3 * sji 

(D22) 

The partial derivatives of the slant 
Z are: 

range, |5^|, 

with respect to, X, Y, and 


1 1 

1 

1 

1 

(D23) 

dx 

- 15^1 


Y - S 2 

(D24) 

ay 

l^i?l 


z - 53 

(D25) 

az 



Thus, the lower row of 


is given by 


3|S^I sISjl 3|S^,I 
ai 9J az 



(Y - S2) 


(Z 



(D26) 


Doppler Radar 

A Doppler radar system is assumed which measures the ground speed, |7|, 
and the angle (drift angle, y) of the ground track from the longitudinal 
axes of the airframe projected in the horizontal plane. These measurements 
are independent of the aircraft position. 


The Jacobian 
defined as: 

measurement matrix, of 

the velocity and drift angle is 



aki 

aki 

aivp 




ax 


3a 



h- 

il 


lx 

(D27) 



_3x 


3s _ 


The groundspeed is 
so that 

the speed of the 

aircraft 

along the reference trajectory 


|V|2 = 

ar^ + H 


(D28) 
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The elements of the upper row of Hp are found by taking the partial 
derivatives of |7l with respect to each of the velocity components. The 
resulting upper row is 



3|F| 

3|7| 

1 

X Y Z 



3a ^ 

r i7i 

L J 


(D29) 


To evaluate the partial derivatives in the lower row of ^Hp, we define 
a new orthogonal coordinate system in which the unit vector U 2 is tangent 
to the ground track, a second unit vector Ui is along the local vertical 
at the aircraft, and a third unit vector is perpendicular to Ui and U 2 , 
such that a right-handed system is^formed.^ The relationship of the ground- 
speed and drift angle to the U\, ^3 coordinate system is shown 

in sketch (Dc) . 



The unit vectors Ui, &2> ^3 found from the following 

relationships : 


A 

Ul 


U3 


h 


R 

CD30) 

1^1 


Fxi? 

CD31) 




(D32) 
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Following the same general arguments as were used to find the partial 
derivatives of the VORTAC/TACAN bearing, we find the lower row of Hr. in the 
^\y U 2 > coordinate system to be: 


9y 9y 9 y 

3^2 9^3^ 



(D33) 


Again following the ardent used with the VORTAC/TACAN bearing, the required 
transformation matrix is the transpose of a matrix which has as its columns 
the unit vectors Ui, U^, and Therefore, the final lower row in the 

% matrix is found using equation (D35) and this transformation matrix as 
follows : 


9y 9y 3y 
_3m 3^ 9s 


= -^[0 1 U2 Uzf 



(D34) 


Radio and Barometric Altimeters 

The radio altimeter and the barometric altimeter (which is compensated 
for in air data computer) both produce altitude and altitude rate information. 

The altitude is given by the magnitude of a vector from the center of 
the earth to the aircraft, R, minus the radius of the earth. R : 

o 

\R\ . R^^ (X^ + y2 22; 1/2 . (.^ 35 ^ 

The altitude rate h is the component of the velocity along the R vector 
that is ' 


• ^ ^ _ XX '+ YY + ZZ 

m m 

The Jacobian matrix is defined as 




dh dh dh 0 0 0 

3Jf 97 92 

3^ 3fe 9fe dfl 3ft 9^ 
97 37 92 37 37 32 


(D36) 


(D37) 
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or 


1 

y 

* 2 
1 I 

r “ “ ” _ ~i 

0 ' 0 ' 0 

(D38) 

V t ^ 

• * Y 

^ ~ ^ \R\ 

% : Z 

\ ^ \R\ \ 

y 1 7 - 2 
' ' 
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APPENDIX E 


FIXED-TIME-OF -ARRIVAL COURSE CORRECTIONS 


Course corrections are intended to return the aircraft to the reference 
path and to arrive there at a specific time. These corrections, which will 
be based on the aircraft estimated state, are computed from a course correc- 
tion law. This course correction law is a fixed -time -of -arrival law which 
allows computation of a commanded change in speed. As, a commanded change in 
aircraft heading, A\p, and a commanded change in aircraft pitch, A0. In this 
implementation, the commanded changes in the first correction maneuver will 
be required to accomplish the desired correction of position error in a fixed 
length of time. At. At the end of the time At, a second correction maneuver 
will be required to aline the aircraft velocity vector along the desired 
reference path. 


Derivation of Av 


Tiie expression for computing the couiiitanded course correction, Av, will 
be derived in this section. This Av is the vector change in velocity which 
will cause the predicted position error to be zero at some time At following 
the execution of the correction. If the correction is the first maneuver in 
the course correction, then At is a fixed interval of time. If the correc- 
tion is the second maneuver, then At is the time-to-go to the destination. 


Assuming the course correction is to take place at time tf,, a transition 
matrix + At, t^) is computed from the equations of motion, the refer- 

ence state, and the estimated altitude rate. Let ^(tj. + At,ti.) be parti- 
tioned in to 3x3 submatrices as follows; ^ ^ 




J>3 $4 


(El) 


Also, let the column vector, x, which is defined at the estimated deviation 
from the reference state, be partitioned into two 3x1 submatrices; 


A , 

x(t^) 




(E2) 


where x\ is the position portion and x^ is the velocity portion. Using 
the upper half of the partitioned transition matrix, the predicted position 
error at + At may be computed from 


86 



Xi(tj^ + Lt) = ■*■ 


^1 

x^(t^y 


CE3) 


where * ^t,t-j^)xi(ty^) and ^z^'^'k ■*■ both predict position 

error at + At. The first matrix product predicts position error due to 

position error at and the second matrix product predicts the position 

error due to velocity error at 

Equation (E3) may be used to find the commanded velocity correction by 
noting that the required velocity increment, when added to X 2 (tj^) , will cause 
the position error to be zero at * At. Thus, equation (E3) becomes 


0 





+ At.y] 


Xi(tj,) 

Jz ''V 


(E4) 


The desired value of the velocity correction is found by solving equation (E4) 
for Ay. The result is 


AV = + At,tk) I I^x(t^) 

= Bx(t^) 


(E5) 


Implementation of the Ay Correction 

Having determined Ay from equation (E5), the equivalent changes in 
heading, Aip, pitch, A0, and the speed correction. As, can be computed. 

To find the angular changes A^ and A6, a coordinate system is defined 
in which the unit ^vector lies along the estimated velocity vector V^, 

the unit vector Uy is perpendicular to the plane of and the radius* 

vector from the center of the earth, and the unit vector is directed 

upward, perpendicular to the plane of and m^, so as to form a right- 

handed coordinate system. This coordinate system is shown in sketch, (Ea) 
which also shows the vector Ay and its components Ay^, Ay^, and Ay^, as 
well as the desired final velocity 7^ in the coordinate system just defined. 
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In this sketch, the angle between Fg and Vf in the u^., Uy plane 
is the change in heading, Aip, and is found from the expression 


tan Alp - 


Ay,. 


\V^\ * Av 


X 


CE6) 


where 


Ay Ay Ay 
X y z 


r 


= Ay 


/s 


U V u 
X y z 


CE7) 


Similarly, the pitch change angle, A6, is found from 


tan A6 = 


Ay; 


|F^| + Ay 


X 


(E8) 


The speed correction (assuming the pitch and heading changes are executed 
first) is given by 


= \v^\ - [V^l = \V^ * AV\ - |F^| 


(E9) 
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Course Correction Execution Error Statistics 

A new coordinate system, > is defined which is the same as 

that shown in the preceding sketch, except that lies along Vj>. This 

coordinate system is shown in sketch (Eb) . In addition, the actual executed 
velocity, is shown along with the random error in measuring the changes in 
the pitch and heading angles, 60 and 6tJ;. 



The vectors Vf and V can be expressed in matrix notation using the 
technique developed*' in reference 11. Thus, after defining 6s to be the 
error in measuring As we have 


= \V^\M 


(ElO) 


and 




"cos 

56 cos 

6ti^ 

7 = ( 

'\VJ + 6s)m 

cos 

66 sin 


e ' 


sin 

66 



(Ell) 


where W is a transformation matrix which relates the coordi- 

nate system to the reference ALT, ALONG-TRACK, CROSS-TRACK coordinate system. 
The angles 56 and 6ip are small angles, so that the following approxima- 
tions are made: 
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cos 66 = 1 


cos 6^ = 1 

sin 66 = 66 

sin 6ip = Sip ^ 


Equation (E12) now becomes 


V 

e 


(\V^\ + Ss}M 


1 

6ip 

66 


The error in making the course correction is 



"l' 


II' 

(\Vj,\ + Ss)M 

Sip 

- \Vf\M 

0 


66j 


.0, 



'o' 



= \Vf\M 

Sip 

+ SsM 



66 




The covariance matrix of the errors in making the correction is found by- 
taking the expected value of nn^. Thus 



0 

0 

0 ^ 


1 

0 

0 

= E(\v^\^)M 

0 

a ^ 
ip 

0 

tf- + cr 

0 

2 

0 


0 

0 

V. 


0 

0 



where the following definitions have been made: 


(E12) 


CE13) 


(E14) 


(E15) 


90 



Er«502; = 
E(^8^) = a ^ 

B 


CE16) 


In addition, it has been assvuned that 60, 6ip, and 6a are independent random 
variables with zero means and that 6s is independent of V^. 

To evaluate E(\Vjy\^}, note that 

I 7^1 2 = (V^ + hv)^(V^ + Ay; (E17) 


After expanding, taking the expected value of both sides, and noting that 7^ 
is a constant vector, we have 

E(\v^\^) = + 2Er7^*Av; + E(hv^^v) (E18) 


The course correction vector. Ay, is assumed to be spherically distributed 
with zero mean. While this assioinption is not particularly good, it does not 
lead to a significant error because of the dominance of the estimated veloc- 
ity vector 7g. The main reason for making the assumption is its simplifica- 
tion of the expressions given below. 

With the above assumption, the expected value of each component of Au 
will have a zero mean. Therefore, 

S’ ('7^- Ay; = 0 CE19) 

To simplify, let the following definition be made: 

fffAy'Ay; = Ay^ (E20) 

rme 

Equation CE19) now becomes 

I 7^1 2; = |7^[2 + AyJ^^ (E21) 

where Ay is the rms computed correction. The quantity Ay can be 

• t * JL*TnS 

found by notxng that 

Ay2 = Trace E(Lv Ay^; (E22) 

rms 

Using equation (E5) , equation (E22) becomes 

Ay2 = Trace BE(xsF (E23) 
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In reference 1, it is shown that 

E(^^) = R(t:j^) - P(tj^) 


so that 

Ay^ = Trace B[R(t-^} - 


(E24) 


(E25) 


where the matrix R(t^) is a covariance matrix of the errors in the devia- 
tions of the actual state from the reference state and is a covari- 

ance matrix of the errors between the actual and estimated states. Using the 
above expressions, equation (E16) may now be written as: 


0 


E(r\rf ) 



+ Trace B 





M 


0 


0 




0 0 



10 0 


+ 


a 

S 








(E26) 


Correction of P(t^) and P(t^) 

Following a course correction, the covariance matrices P and R must 
be corrected for the random errors which were generated by the course correc- 
tion process. TTie basic expressions for this correction are derived in 
reference 1 and will not be repeated here. Instead, the final results will 
be discussed. 


The expression for correcting the PCtz) matrix for the added errors 
resulting from a course correction is : ^ 


p(t^) = p(tp . 


0 srnri ; 


(E27) 


where all submatrices are 3x3* 



To correct R(f^), first let 


G = 






(E28) 


where B 
by: 


is 3x6 and G is 6x6, Then the expression for 
R(t^) = r j + G)^(t^) - + c?;^ + p(t^) 


is given 

(E29) 
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